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Affine  invariant  interdirection  tests  are  proposed  for  a repeated  measures  problem. 
The  test  statistics  proposed  are  applications  of  the  one-sample  interdirection  sign  test  and 
interdirection  signed-rank  test  to  a repeated  measurement  setting.  The  interdirection  sign 
test  has  a small  sample  distribution-free  property  and  includes  the  two-sided  univariate  sign 
test  and  Blumen’s  bivariate  sign  test  as  special  cases.  The  interdirection  signed-rank  test 
includes  the  two-sided  univariate  Wilcoxon  signed-rank  test  as  a special  case.  The 
asymptotic  null  distributions  of  the  proposed  statistics  are  obtained  for  the  class  of 
elliptically  symmetric  distributions.  In  addition,  the  asymptotic  distributions  of  the 
proposed  statistics  under  certain  contiguous  alternatives  are  obtained  for  elliptically 
symmetric  distributions  with  a certain  density  function  form.  Comparisons  are  made 
between  the  proposed  statistics  and  several  competitors  via  Pitman  asymptotic  relative 
efficiencies  and  Monte  Carlo  studies.  The  interdirection  tests  proposed  appear  to  be  robust. 
The  sign  test  performs  better  than  the  other  competitors  when  the  underlying  distribution  is 
heavy-tailed  or  skewed.  For  normal  to  light-tailed  distributions,  the  Hotelling’s  T2  and 
signed-rank  test  have  good  powers  when  the  variance-covariance  structure  of  the 


IX 


underlying  distribution  is  non  H-type,  otherwise  ANOVA  F and  the  rank  transformation 
test  RT  perform  better  than  the  others. 

An  alternative  test  for  the  one-sample  multivariate  location  problem  is  also  proposed 
which  extends  the  univariate  signed-rank  test  to  multivariate  settings.  The  test  proposed  is 
somewhat  like  applying  the  interdirection  sign  test  to  the  sums  of  pairs  of  observed  vectors. 
It  includes  the  two-sided  univariate  Wilcoxon  signed-rank  test  as  a special  case.  The 
asymptotic  distributions  of  the  proposed  statistic  under  the  null  hypothesis  and  under 
certain  contiguous  alternatives  are  obtained  for  a class  of  elliptically  symmetric 
distributions.  Comparisons  are  made  between  the  proposed  statistic  and  Hotelling’s  T2  via 
Pitman  asymptotic  relative  efficiencies.  The  signed  sum  test  proposed  performs  better  than 
Hotelling’s  T2  when  the  underlying  distribution  is  heavy-tailed.  However,  for  normal  to 
light-tailed  distributions,  the  Hotelling’s  T2  performs  slightly  better  than  the  proposed  test 
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CHAPTER  1 
INTRODUCTION 


In  this  dissertation  we  investigate  test  statistics  for  certain  repeated-measures  and 
one-sample  multivariate  location  problems.  For  the  repeated-measures  problem,  we  let 
X!,  ....  Yn  be  independently  and  identically  distributed  as  Y , where  Y = (Y^  Yp)T is 

from  a p-dimensional,  p > 2,  absolutely  continuous  population.  For  each  subject  i,  we  shall 
regard  Yj  as  repeated  measures  with  one  observation  for  each  of  the  p treatments.  We  use 

the  general  mixed  model  with  a subject  by  treatment  interaction.  In  vector  form,  we 
consider  the  model  (see,  e.g.,  Winer  [1971],  p.  278), 

Yi  = pilp+i  + fiii+£i,i=  1, ...,  n,  (1.1.1) 

where  Ip  is  the  pxl  vector  of  l’s,  1 = (Xj, ...,  xp)T  is  the  vector  of  fixed  treatment  effects, 
and  random  variables  (3j,  £lj,  £j  for  i = 1, ...,  n are  all  mutually  independent.  (More  details 

about  this  model  will  be  given  in  section  2.1.)  Note  that  the  variance-covariance  matrix  of 
fili  is  probably  not  H-type,  to  be  described  later.  We  are  concerned  with  the  test  of  equal 

treatment  effects,  which  in  model  (2.1.1)  can  be  described  as: 

H0  : Tj  =x2=  •••  =Tp  versus  Ha:  Xj  * Xj'  for  some  j * ].  (1.1.2) 


We  first  consider  several  parametric  statistics  for  this  problem  when  the  population 
is  p-variate  normal.  Probably  the  most  well-known  parametric  procedure  for  this  problem 
is  based  on  Hotelling  (1931)  T^  test.  Define  the  (p-l)-variate  random  vector  Zj  by 


— i (Zjj,  ....  Zj  p.j)  - (Yj!  Yjp,  ....  Yjp.2-Yip)T,  i = 1, ....  n. 


1 


2 


The  test  of  (1.1.2)  can  be  carried  out  by  the  Hotelling’s  T2  statistic  computed  from  the 
mean  vector  and  sample  variance-covariance  matrix  of  the  vectors  of  differences  Zj’s.  This 

is  due  to  Hsu  (1938).  The  Hotelling’s  T2  is  then  defined  as 

T2  = T2(  Zh  ...,  Zn  ) = n Z7 (I)'1!,  (1.1.3) 


where 

2 Z Zi,  and  £ = -py  £ (Zi-Z)^-!)7. 
i=l  11  i=l 

Under  Hq,  T is  asymptotically  Chi-square  with  p-1  degrees  of  freedom.  If  the  underlying 
population  is  p-variate  normal,  then  the  null  distribution  of  T2  is  a multiple  of  the 
F-distribution  with  p-1  numerator  degrees  of  freedom  and  n-p+1  denominator  degrees  of 

r\ 

freedom.  Hotelling’s  T is  invariant  with  respect  to  nonsingular  linear  transformations  of 
the  observations  Zj,  i = 1, ...,  n.  That  is,  if  D is  any  nonsingular  (p-l)x(p-l)  matrix,  then 

T2(DZ1,...,DZn)=T2(Z1,...,Zn).  (1.1.4) 


We  shall  call  this  invariance  property  affine-invariance.  This  appealing  invariance  property 
ensures  that  the  value  of  the  test  statistic  remains  unchanged  following  rotations  of  the 
observations  about  the  origin,  reflections  of  the  observations  about  a (p-2)-dimensional 
hyperplane  through  the  origin,  or  changes  in  scale.  Hence  the  performance  of  Hotelling’s 
T test  will  not  depend  on  the  structure  of  the  population  variance-covariance  matrix  or  the 
direction  of  shift. 

Another  parametric  procedure  for  this  problem  is  the  classic  ANOVA  F test.  The 
test  is  based  on  the  original  observations  Yjj,  1 < i < n,  1 < j < p,  and  is  defined  as 

P _ _ 


F = 


n I (Y.j-Y  y 
i=l 


iTT  i,  £ (%Yi.-Y.j+Y. ,)2 

1=1  J=1 


(1.1.5) 
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where  Y ; = - £ Y; 

n i=l 


•J 


U’ 


IP-  - IP 

= - Z Y and  Yj  = - £ Yj:.  If  the  underlying  population  is 
P j=l  'J  ‘ P j=i  J 


p-variate  normal  with  variance-covariance  matrix  of  the  form 


1 = 


(1.1.6) 


then  the  null  distribution  of  F is  F-distribution  with  p- 1 numerator  degrees  of  freedom  and 
(n-l)(p-l)  denominator  degrees  of  freedom.  This  test  for  (1.1.2)  under  the  variance- 
covariance  matrix  (1.1.6)  was  obtained  by  Wilks  (1946)  from  the  generalized  likelihood 
ratio  principle  when  the  underlying  population  is  p-variate  normal.  The  matrix  in  (1.1.6)  is 
said  to  have  compound  symmetry.  While  compound  symmetry  is  a sufficient  condition  for 
the  test  statistic  F to  have  an  exact  F-distribution,  it  is  not  a necessary  one.  Huynh  and 
Feldt  (1970)  have  found  a necessary  and  sufficient  condition,  which  may  be  expressed  in 
three  alternative  forms  (see  Morrison  [1976],  p.152): 

(1)  The  population  variance-covariance  matrix  L = (a-,),  1 < j,j'  < p,  has  the  pattern 


defined  by 


JJ 


aj+aj-U  if  j = j' 
OCj-HXj.  if  j * j' 


(1.1.7) 


where  X,  ah  ...,  ap  are  (p+1)  arbitrary  constants  such  that  the  resultant  matrix  is  positive 
definite. 

(2)  All  possible  differences  Yj-Yj1  (j  ^ j')  of  the  response  variates  have  the  same  variance. 

(3)  Define  e,  a function  of  the  elements  of  Z,  by 
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P2(°d  - <L.)2 


(1.1.8) 


e = 


e = 


(P-I)(v  V + p2(j2  ) 


Then  an  alternative  statement  of  the  necessary  and  sufficient  condition  is  that  e = 1. 

It  may  be  noted  that  the  scalar  factor  e is  the  multiplicative  adjustment  factor  for  the  degrees 
of  freedom  proposed  by  Box  (1954)  and  by  Greenhouse  and  Geisser  (1959).  Any  matrix 
whose  elements  satisfy  (1.1.7)  is  referred  to  as  a matrix  of  Type  H.  If  the  underlying 
population  is  p-variate  normal  with  H-type  variance-covariance  matrix,  the  power  of  the  F 
test  will  always  be  greater  than  that  of  the  Hotelling’s  T2  test  for  the  same  alternative 
because  the  second  degrees-of-freedom  parameter  of  the  F test  is  greater  than  that  of 
Hotelling’s  T2  test.  However,  unlike  Hotelling’s  T2,  ANOVA  F does  not  have  affine- 
invariance.  Thus,  the  performance  of  F test  depends  strongly  on  the  structure  of  the 
variance-covariance  matrix  of  the  underlying  population. 

Morrison  (1972)  proposed  several  tests  for  this  problem  under  various  assumptions 
on  the  variance-covariance  matrix.  The  tests  include  the  usual  Hotelling-Hsu  statistic  and 
the  classical  ANOVA  F test  as  special  cases. 

Many  nonparametric  competitors  to  Hotelling’s  T2  have  been  proposed.  A well- 
known  statistic  is  a rank  test  due  to  Friedman  (1937).  His  test  statistic  is  based  entirely  on 
within-block  information  and  ignores  between-block  information.  Let  Rjj  denote  the  rank 
assigned  to  Yy  within  subject  i.  Let  Rj  represent  the  sum  of  the  ranks  associated  with 


treatment  j,  i.e.,  Rj  = £ Rjj,  j = 1,  ...,  p.  The  Friedman  test  statistic  for  this  case  with  no 


n 


(1.1.9) 
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Under  H0,  the  test  S is  asymptotically  Chi-square  with  p-1  degrees  of  freedom.  (See,  e.g., 
Hollander  and  Wolfe  [1973],  p.140.)  Define 

(n-l)S 
_n(p-l)-S  • 

Then  under  H0,  the  test  Fs  is  compared  with  an  F-distribution  with  degrees  of  freeedom 
p-1  and  (n-l)(p-l).  A more  accurate  approximation  to  Friedman’s  test  is  proposed  by 
Jensen  (1977).  Iman  and  Davenport  (1980)  proposed  two  new  approximations  and  also 
pointed  out  that  the  F approximation  is  better  than  the  %2  approximation. 

Another  rank  test  was  proposed  by  Koch  (1969).  His  test  statistic  used  the  ranks  of 
the  aligned  observations,  obtained  by  subtracting  from  each  observation  the  average  of  the 
observations  for  that  block.  This  alignment  process  will  eliminate  or  at  least  reduce  the 
block  effect.  To  introduce  his  test,  let  us  define 

Rjj  = rank  of  (Yy-Yj)  among  Yn-YL, ...,  Ynp-Yn-, 


R.J  = i|1RiiandRi.  = f.|Rij- 


Then  the  test  statistic  has  the  form 


nf  (Rj  - n-^1)2 


W‘  = EL_1 i 


1 n p 


n(p-l)  £ £ ^1J 


d \2 


I X (Rij  - Ri.) 


(1.1.10) 


If  Yj,  ...,  Yn  are  i.i.d.,  then  under  the  null  hypothesis  of  the  exchangeability  of  the 
components  of  Yj  = (Yjj,  ...,Yjp)T,  the  test  W has  an  asymptotic  Chi-square  distribution 
with  degrees  of  freedom  p- 1 . 

Unlike  the  Friedman  test,  which  depends  entirely  on  within-block  rankings,  Quade 
(1979)  considered  a p-sample  extension  of  the  Wilcoxon  signed  rank  test  by  taking 


6 


advantage  of  the  between-block  information.  This  is  done  by  considering  weights  assigned 
to  each  block  on  the  basis  of  some  measure  of  within-block  sample  variation,  such  as  the 
range,  standard  deviation,  mean  deviation,  or  interquartile  difference.  To  illustrate  this  test 
statistic,  let  us  define 

Dj  = D(Yj),  a location-free  statistic  that  measures  the  variability  within  the  ith  block, 
Qi  = rank  of  D,  among  Db  . . .,  Dn,  and 


Ry=rank  of  Yy  among  Yn, ...,  Yip. 


Quade's  procedure,  based  on  weighted  within-block  ranks,  is  defined  as 


W 


P r n -1  r\ 

72l[lQiRij]2 

j=l  i=l 

p(p+l)n(n+l)(2n+l) 


9(p+l)n(n+l) 

2(2n+l) 


72  i [ £ QifRij-^)]2 

1=1  i=l 

p(p+l)n(n+l)(2n+l)  ‘ 


(1.1.11) 


Under  the  null  hypothesis  of  the  exchangeability  of  the  components  of  Yj,  the  test  W has  an 
asymptotic  Chi-square  distribution  with  p-1  degrees  of  freedom. 

A natural  nonparametric  analog  of  ANOVA  F test  was  proposed  by  Iman,  Hora  and 
Conover  (1984).  Their  procedure  is  first  to  transform  all  the  observations  to  ranks  from  1 
to  np  and  then  apply  the  parametric  ANOVA  F test  to  the  ranks.  This  approach  retains  both 
the  within-  and  between-block  information.  Defining 

Rjj  = rank  of  Yy  among  Yn, ...,  Ynp, 


the  rank  transformation  test  is  defined  as 
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RT  = 


nf  (R.-R  Y 
_ >1 


n-1  .? ;(Rij  Ri-  R.j+R..)2 

i=l  j=l 


(1.1.12) 


1 " 


1 P 


1 P 


where  R j - n X Rij,  Ri.  = ~ Z Rij  and  R = - Z R.j-  Assuming  Yy's  are  mutually 
independent,  and  considering 

H0  : Fij  = Fifor  i = 1, ....  nandj  = 1,  ...,  p, 

where  Fy  is  the  distribution  function  of  Yy,  then  the  asymptotic  null  distribution  of 
(p-l)RT  is  Chi-square  with  p-1  degrees  of  freedom  under  suitable  conditions.  Their 


simulations  showed  that  the  behavior  of  test  RT  is  closely  approximated  by  the 
F-distribution  with  (p-1)  and  (n-l)(p-l)  degrees  of  freedom.  Comparisons  made  among 
ANOVA  F test,  Friedman’s  test  (S),  Quade's  test  (W),  and  rank  transformation  test  (RT) 
via  Monte  Carlo  studies  showed  that  the  F test  had  the  most  power  for  normal 
distributions,  the  Quade's  test  and  F test  were  almost  equivalent  and  gave  the  best  results 
for  uniform  distribution,  the  Friedman's  test  and  the  RT  test  gave  similar  results  and  were 
best  for  the  Cauchy  distribution,  and  the  RT  test  has  the  most  power  for  double  exponential 
and  lognormal  distributions.  Hora  and  Iman  (1988)  developed  the  limiting  noncentrality 
parameters  of  the  rank  transformation  statistic  and  some  other  tests,  which  were  then 
evaluated  to  make  comparisons  among  those  tests  via  Pitman  asymptotic  relative 
efficiencies. 

Agresti  and  Pendergast  (1986)  also  considered  a test  that  is  appropriate  when  the 
null  hypothesis  (1.1.2)  is  expressed  as  the  exchangeability  of  the  components  of  Yp  Their 

procedure  utilizes  a single  ranking  of  the  entire  sample.  Let 


Rij  = rank  of  Yy  among  Yn, ...,  Ynp,  Rj  = i £ Rijt 
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p = Corr(Rjj,  Rjj-),  j * j',  and  ct2  = Var(Rjj). 


The  test  statistic  is  based  on 


(1.1.13) 


o (1-P) 


This  test  includes  Koch's  test  and  the  rank  transformation  test  as  special  cases.  They 
argued  that  under  Hq  their  statistic  has  an  asymptotic  Chi-square  distribution  with  p-1 
degrees  of  freedom  if  the  asymptotic  distribution  of  R = (R  A,  ...,  R p^)7  is  (p-l)-variate 
normal.  Since  they  did  not  present  conditions  guaranteeing  this  normality,  Kepner  and 
Robinson  (1988)  concluded  the  work  of  Agresti  and  Pendergast  by  determining  reasonably 
sufficient  conditions  for  R to  have  a (p-l)-variate  normal  limiting  distribution. 

In  analogy  to  the  Iman,  Hora  and  Conover  (1984)  proposal  of  a rank  transformed 
version  of  ANOVA  F test,  Agresti  and  Pendergast  (1986)  considered  a rank  transformed 
version  of  Hotelling's  test.  This  procedure  is  appropriate  when  the  hypothesis  of  no 
treatment  effects  is  more  broadly  expressed  as  the  marginal  homogeneity  condition 
F1  = f2  = •••  = Fp>  where  Fj, ...,  Fp  denote  the  one-dimensional  marginal  distribution  of 

*T* 

Y = ( Y j , . . . , Yp)  . Their  statistic  is  based  on 


RTh  = n RTS_1R, 


(1.1.14) 


where 


i=l 


Their  simulations  showed  treating  the  null  distribution  of  RTH  as  a multiple  of  F- 

distribution  with  p-1  numerator  degrees  of  freedom  and  n-p+1  denominator  degrees  of 
freedom  is  a reasonable  approximation.  They  also  argued  that  under  H0,  their  statistic  has 
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an  asymptotically  Chi-square  distribution  with  p- 1 degrees  of  freedom  if  the  asymptotic 
distribution  of  (R  j,  Rp_i)T  is  (p-l)-variate  normal.  Their  simulations  showed  that  the 
RT  and  RTH  statistics  behaved  much  like  their  parametric  analogs. 

We  now  consider  the  multivariate  tests  for  the  one-sample  location  problem.  As  a 
first  step,  we  let  Xj,  ...,  Xn  be  i.i.d.  as  X = (Xj,  ...,  Xp)T,  where  X is  from  a p-variate 

absolutely  continuous  population  with  location  parameter  §*  (pxl).  We  would  like  to  test 

Ho:£*  = 0 versus  Ha:£**0.  (1.1.15) 

Here  0 is  used  without  loss  of  generality,  since  H0  : 0*  = 0Q  can  be  tested  by  subtracting 
00  from  each  observation  X;  and  testing  whether  these  differences  (X;  - GqI’s  are  located  at 

Q. 

The  classical  procedure  used  in  this  setting  is  Hotelling’s  T2,  which  is  defined  as 
T2  = T2(  Xh  . . .,  Xn  ) = n X1©-^, 

where 

X=-^-ixi;  and£=-prI(Xi-X)(Xi-X)T. 

If  the  underlying  population  is  p-variate  normal,  then  the  null  distribution  of  T2  is  a 
multiple  of  the  F-distribution  with  p numerator  degrees  of  freedom  and  n-p  denominator 
degrees  of  freedom.  The  affine-invariant  property  of  Hotelling’s  T2  is  discussed  earlier  on 
page  2. 

Many  nonparametric  procedures  have  been  proposed.  The  most  popular  statistic  is 
the  component  sign  test,  which  is  a nonparametric  analog  of  Hotelling’s  T2  using  signs  of 
Xjj’s,  1 < i < n,  1 < j < p.  (See,  e.g.,  Randles  [1989],  p.  1045.)  Let 

Sj  = X sgn(Xij),  and£=  (Sb  ...,  Sp)T, 
i=l 
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where  sgn(x)  = 1(0,  -1)  for  x >(=,  <)  0.  The  test  is  based  on 

S*  = S.T(  n$f  )'1S,  (1.1.16) 

where  $/=  (wjj-),  wjj-  =-jpI  sgn(Xij)sgn(XijO,  for  1 <j,j'  <p. 

$ 

Under  Hq,  the  test  statistic  Sn  is  asymptotically  Chi-square  with  p degrees  of  freedom.  Sen 

and  Puri  (1967)  considered  other  score  function  versions  of  this  statistic  as  well.  Bennett 
(1962)  considered  a similar  sign  test  for  the  problem  of  testing  the  equality  of  location 
parameters  in  two  p-variate  (p  < 4)  populations.  Chatteijee  (1966)  also  studied  a similar 
procedure  for  a bivariate  case.  In  the  paper  of  Bickel  (1965),  he  proposed  several  general 
test  statistics  based  on  Hotelling’s  T2.  They  are  Hotelling’s  tests  of  type  I,  M 2 and  W2, 
and  Hotelling’s  tests  of  type  II,  if  2 and  $2.  His  p-variate  test  statistics  M 2 and  fo  2, 
and  W 2 and  t^2  are  quadratic  forms  involving  coordinate-wise  sign  and  Wilcoxon  signed 
rank  statistics,  respectively. 

Unlike  Hotelling’s  T2,  the  above  competitors  to  T2  are  not  invariant  under 
nonsingular  linear  transformations  of  the  observations.  Consequently,  their  power  and 
efficiency  depend  on  the  direction  of  shift  and  the  variance-covariance  matrix  of  the 
alternative  distribution.  In  an  effort  to  overcome  this  problem,  Dietz  (1982)  proposed  sign 
and  signed  rank  tests  for  bivariate  location  settings.  Her  tests  are  three-step  procedures. 
First,  a certain  transformation  is  applied  to  the  observations.  Second,  new  coordinate  axes 
are  chosen.  Third,  standard  sign  and  signed  rank  tests  are  performed  using  the  transformed 
and  rotated  observations.  In  Dietz  (1984),  she  extended  those  results  to  general  linear 
signed  rank  tests  for  multivariate  location  problems.  Here,  however,  a two-stage  procedure 
is  recommended,  corresponding  to  steps  one  and  three  of  the  bivariate  case.  The  resulting 
test  statistics  are  no  longer  invariant  under  linear  transformations,  but  for  elliptically 
symmetric  alternatives,  their  asymptotic  efficiencies  are  independent  of  the  direction  of  shift 
and  the  variance-covariance  matrix.  Her  simulation  study  showed  that  linear 
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transformations  have  little  effect  on  the  small  sample  power  of  the  tests  for  nearly 
degenerate  distributions. 

Two  well-known  bivariate  sign  tests  are  due  to  Hodges  (1955)  and  Blumen  (1958). 
Their  procedures  are  affine-invariant  and  have  a distribution-free  property.  Joffe  and  Klotz 
(1962)  presented  an  expression  for  the  exact  null  distribution  of  the  Hodges  bivariate  sign 
test.  They  also  computed  the  Bahadur  limiting  efficiency  of  the  test  relative  to  the 
Hotelling’s  T2  test  for  normal  alternatives.  Killeen  and  Hettmansperger  (1972)  made  an 
exact  Bahadur  efficiency  comparison  of  Hotelling’s  T2  with  respect  to  both  Hodges’  and 
Blumen's  bivariate  sign  tests.  Klotz  (1964)  obtained  exact  power  for  the  bivariate  sign  tests 
of  Hodges  and  Blumen  under  normal  alternatives  and  therefore  permitted  comparisons  of 
the  two  tests  for  sample  sizes  n = 8 through  12. 

The  procedure  proposed  by  Bennett  (1964)  for  the  bivariate  case  is  a signed-rank 
test  generalizing  Wilcoxon's  univariate  signed-rank  test.  This  test  is  not  affine-invariant. 

Another  affine-invariant  bivariate  rank  test  was  introduced  by  Brown  and 
Hettmansperger  (1985).  Their  statistic  is  based  on  the  gradient  of  Oja's  measure  of  scatter 
(Oja,  1983).  Letting  A(Xj,  Xj;  £ ) denote  the  area  of  the  triangle  formed  with  Xj,  Xj  and 
£*  as  vertices,  define 


T(fi*)  = ZIA(Xi,Xj;fi*).  (1.1.17) 

i<j 

This  is  the  Oja  measure  of  scatter.  The  value  which  minimizes  T(0*)  is  the  Oja 
generalized  median  of  the  bivariate  sample.  Brown  and  Hettmansperger  proposed  the  use 
°f  Qn(fi  ) = dT((i  ),  the  vector  of  partial  derivatives  of  T(0  ).  The  generalized  median  of 
Oja  is  the  value  which  minimizes  I Qn(£*)  I.  Define 

S = £ Q2n(-Xi),  and  C = £ Q2n(-Xi)Q2n(Xi), 


(1.1.18) 
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where  Q2n(-Xi)  is  computed  using  observations  Xj, Xn,  and  their  projections  through 

the  origin,  -Xj,  -Xn.  Their  test  for  (1.1.15)  is  defined  as  S7^'1^.  Under  Hq,  the  test 

statistic  is  asymptotically  Chi-square  with  2 degrees  of  freedom. 

Oja  and  Nyblom  (1989)  also  studied  the  bivariate  location  problem.  Their  tests  are 
analogs  of  the  univariate  sign  test.  Denote  the  direction  angle  of  Xj  by  9,.  Then  6*  = 0J+7U 

/ r / 

(mod  27t)  is  the  direction  angle  of  -Xj.  Write  < 02  < ...  < 02n  for  the  ordered  angles  in 
the  set  {0j,  ...,  0n,  0 x , ....  0n).  Define 

Zj  = 1 ifO'e  {0j 0n) 

0 if0je  {0*,  ...,  0*},i=l,  ...,  n,  (1.1.19) 

and 

Zn+i  = 1-Zj , i =1,2,....  (1.1.20) 

T* 

The  vector  Z = (Zj,  ...,7^)  indicates  which  of  the  observations  lie  above  or  below  the 
horizontal  axis.  They  proposed  using  test  statistics  based  on  Z.  The  test  statistics  are 
distribution-free  and  affine-invariant,  and  include  Hodges  (1955)  and  Blumen  (1958)  sign 
tests  as  special  cases.  Also,  they  proposed  some  new  intuitively  appealing  tests.  A general 
class  of  these  invariant  sign  test  statistics  is 

Z [£  (Zk+i-Hh(i/n)]2  (1.1.21) 

k=0  i=l  z 

where  h is  a suitably  chosen  score  function. 

In  the  recent  study  Randles  (1989)  proposed  an  interdirection  sign  test  for  this 
problem.  His  test  statistic  included  the  two-sided  univariate  sign  test  and  the  Blumen 
(1958)  bivariate  sign  test  as  special  cases.  Also,  Peters  and  Randles  (1990)  suggested  a 
signed-rank  test  based  on  interdirections,  which  includes  the  two-sided  Wilcoxon  signed- 
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rank  test  as  a special  case.  The  bivariate  case  of  their  statistic  was  considered  in  the 
dissertation  of  Peters  (1988).  The  interdirection  sign  test  and  the  interdirection  signed-rank 
test  will  be  described  in  detail  in  Chapters  2 and  3,  respectively,  where  they  are  applied  to  a 
repeated  measures  problem. 

In  this  dissertation,  the  interdirection  sign  test  for  a repeated  measures  problem  is 
defined  in  Chapter  2.  The  asymptotic  distributions  of  the  test  under  H0  and  under  certain 

contiguous  alternatives  are  obtained  in  sections  2.2  and  2.3,  respectively.  The  Pitman 
asymptotic  relative  efficiencies  of  the  test  relative  to  Hotelling’s  T2  are  presented  in  the  last 
section.  In  Chapter  3,  the  interdirection  signed-rank  test  for  the  same  problem  is  described, 
with  its  asymptotic  distributions  obtained  in  section  3.2  and  the  evaluations  of  the  ARE  of 
the  signed-rank  test  relative  to  Hotelling’s  T2  established  in  the  section  3.3.  Comparisons 
of  several  competing  procedures  are  made  in  Chapter  4 via  Monte  Carlo  studies.  An 
alternative  test  for  the  one-sample  multivariate  location  problem  is  proposed  in  Chapter  5. 
Some  useful  intermediate  results  are  presented  in  section  5.2.  The  asymptotic  distributions 
of  the  test  under  Hq  and  under  certain  contiguous  alternatives  are  developed  in  sections  5.3 

and  5.4,  respectively.  Finally,  in  section  5.5,  we  evaluate  the  ARE  of  the  proposed  test 
relative  to  Hotelling’s  T2. 


CHAPTER  2 

A MULTIVARIATE  SIGN  TEST  BASED  ON  INTERDIRECTIONS 
FOR  REPEATED-MEASURES  DESIGNS 


2, 1 Definition  of  the  Test  Statistic 

The  multivariate  sign  test  based  on  interdirections,  denoted  by  Vn,  was  proposed  by 
Randles  for  the  one-sample  multivariate  location  problem.  In  this  section  we  will  show 
how  this  test  statistic  can  also  be  applied  to  repeated-measures  designs  for  detecting 
treatment  effects. 

For  the  one-sample  multivariate  location  problem,  we  let  Xlf  ...,  Xn,  where  X{  = 

(Xji,  •••,  Xip)T,  be  independent  and  identically  distributed  (i.i.d.)  as  X = (X^  ...,  Xp)T, 

where  X is  from  a p-dimensional  absolutely  continuous  population  with  location  parameter 
$ 

Q (pxl).  We  would  like  to  test 

Ho  : 0*  = 0 versus  Ha:  0*  * 0. 

Here  0 is  used  without  loss  of  generality,  since  Hq  : 9 = 0q  can  be  tested  by  subtracting 
00  from  each  observation  Xj  and  testing  whether  these  differences  (Xj  - 0q)’s  are  located  at 
Q. 

For  the  problem  of  single-factor  repeated-measures  designs,  we  let  Yj, ...,  Yn  be 
i.i.d.  as  Y = (Yj, ...,  Yp)^,  where  Y is  from  a p-dimensional,  p > 2,  absolutely  continuous 
population.  Note  that  the  components  of  Yj  are  repeated  measurements  of  the  ith 

experimental  unit.  We  will  use  the  general  mixed  model  with  a subject  by  treatment 
interaction.  In  vector  forms,  we  consider  the  model  (see,  e.g.,  Winer  [1971],  p.  278), 

Xj  = Pilp  + 1 + fili  + £j,  i = 1,  ...,  n,  (2.1.1) 
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where  lp  is  the  pxl  vector  of  l’s,  l = (ij, xp)T  is  the  vector  of  fixed  treatment  effects. 
Pi  represents  the  random  effect  of  the  ith  subject,  Pij  denotes  the  vector  of  the  ith  subject 
by  treatments  interactions,  and  e,  is  the  vector  of  random  error  of  the  ith  subject.  We 
assume  pj’s  are  i.i.d.  with  mean  0,  j3i[’s  are  i.i.d.  with  mean  0 and  a general  variance- 

covariance  matrix  possibly  not  H-type,  described  in  Chapter  1,  gj’s  are  i.i.d.  with  mean  0 
and  variance-covariance  matrix  aeIp,  where  Ip  is  the  pxp  identity  matrix.  The  random 

variables  pj,  Pi,,  £,  for  i = 1,  ...,  n are  all  mutually  independent.  We  are  concerned  with 
the  test  of  equal  treatment  effects,  which  in  model  (2.1.1)  can  be  described  as: 

Hq  : Xi  = X2  = ...  = xp  versus  Ha:  Xj  * Xj'  for  some  j j . (2.1.2) 


Note  that  the  sample  Yj,...,  Yn  has  location  parameter  i.  The  problem  is  to  test  whether 
the  components  of  the  location  parameter  j are  all  equal.  We  can  transform  this  problem  to 
the  standard  one-sample  multivariate  location  problem,  described  earlier,  by  looking  at  the 
differences  among  the  components  Y;j  within  each  observation  Yj.  The  transformation  is 
described  below.  Define 


Yjl 

Yi2-Yip 


-1 

-1 

-i 


Yii 

Yi2 


Yj  p-i 

- Yip  - 


= A Yj , i = 1, ...,  n, 


and 


" Xi-tp  " 

*2-^p 

■ i -r 

x2 

• 

= 

i -i 

’ i -i 

• 

-Vi“V 

VI 

-^p J 
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where 


1 


A = 


-1 

-1 

-i 


(2.1.3) 


is  a (p-l)xp  matrix.  Now,  we  have  a (p-l)-variate  sample  Zt,  ....  Z^,  which  can  be 
modeled  via 

Z;  = 0 + +8*  , i = 1,  ...,  n. 


where  0 ((p-l)xl)  is  the  location  parameter.  Note  that  the  variance-covariance  matrix  of 
* , 

pij  is  possibly  not  H-type.  Thus,  testing  the  hypotheses  given  in  (2.1.2)  is  equivalent  to 
testing 

Hq:S  = Q versus  Ha:0*O.  (2.1.4) 

We  have  shown  that  the  test  of  (2.1.2)  based  on  p-variate  observations  Yj, ...,  Yn  can  be 

earned  out  by  using  a multivariate  location  test  based  on  a statistic  like  the  interdirection 
sign  statistic  Vn,  computed  on  the  transformed  (p-l)-variate  observations  Zh  ...,  Z for 
testing  (2.1.4).  Now  we  will  describe  the  test  that  rejects  H0  : 0 = 0 for  large  values  of  the 
statistic 

Vn  I Icosfrp^),  (2.1.5) 

i=l  k=l 


where 


ifi*k 


0 


if  i = k, 


(2.1.6) 
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- *[(£)-($]• 

and  Cj^  denotes  the  number  of  hyperplanes  formed  by  the  origin  0 and  other  p-2 
observations  (excluding  Z,  and  Zk)  such  that  Zj  and  Z^  are  on  opposite  sides  of  the 
hyperplane  formed.  The  counts  {C^  I 1 < i < k < n},  called  interdirections,  are  used  via 
7ipik  to  measure  the  angular  distance  between  Zj  and  Zk  relative  to  the  positions  of  the 
other  observations.  This  statistic,  Vn,  includes  Blumen’s  bivariate  sign  test  and  the  2-sided 

univariate  sign  test  as  special  cases.  Also,  it  is  affine-invariant  and  has  a distribution-free 
property  under  Hq,  for  a broad  class  of  distributions,  called  distributions  with  elliptical 

directions,  which  includes  all  elliptically  symmetric  populations  and  many  skewed 
populations  as  well.  In  the  next  section,  we  will  concentrate  on  the  family  of  elliptically 
symmetric  populations. 


2.2  Null  Distribution  of  Vn 

In  this  section  we  will  find  the  null  distribution  of  Vn  under  the  class  of  elliptically 
symmetric  distributions,  which  is  defined  below. 

Definition  2.2.1  Assuming  the  existence  of  a density  function,  the  mxl  random  vector  X 
is  said  to  have  an  elliptically  symmetric  distribution  with  parameters  ji  (mxl)  and  Z (mxm) 
if  its  density  function  is  of  the  form 

fx(x)  = KmIZr1/2h[(x  - fctiVfx  - U)],  (2.2.2) 

for  some  non-negative  real-valued  function  h,  where  Z is  positive  definite  and  Km  is  a 
scalar.  We  will  write  this  distribution  of  X as  Em(ji,  Z). 
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Throughout  this  chapter,  we  will  use  Yj, Yn  to  denote  the  original  sample,  and 
use  Zj,  Zn,  defined  in  (2.1.3),  to  denote  the  transformed  sample,  to  which  the  test 
statistic  Vn  is  applied.  Now,  let’s  assume  Yj,  ...,  Yn  are  i.i.d.  as  Y = (Yj,  Yp)T, 
where  Y is  Ep(ji,  Z).  To  apply  the  result  of  null  distribution  of  Vn  under  elliptically 

symmetric  distributions,  proved  by  Randles,  we  shall  first  prove  that  the  transformed 
sample  Z\ , ....  Zn  is  also  elliptically  symmetric.  To  do  this,  we  use  the  following  lemma, 

which  was  given  as  an  exercise  in  Muirhead’s  (1982)  book. 

Lemma  2.2.3  If  X is  Em(y.,  Z)  then  : 

•jfTy 

(i)  the  characteristic  function  <J>x(t)  = E(e  - — ) has  the  form 

• T 

$x(t)  = e1-  ^ y(tTZt)  for  some  function  \j /,  (2.2.4) 

and 

(ii)  provided  they  exist,  E(X)  = g and  Cov(Y)  = aS  for  some  constant  a. 

Theorem  2.2,5  If  Y_  is  Ep(ii,  Zj  and  Z = A Y,  defined  in  (2.1.3),  then  Z is 
Ep.i(A  y.,  A Z At). 

Proof  of  Theorem  2.2.5  Since  Y is  Ep(ji,  Z),  by  Lemma  2.2.3,  the  characteristic  function 
of  Y at  I,  a pxl  vector,  has  the  form  of 

• T 

<t>x(i)  = ^ V(lTS  t)  for  some  function  \\ f. 

Thus,  the  characteristic  function  of  Z at  s,  a (p-l)xl  vector,  is 
bz(s)  = E(ei-T— ) = E(ei-T— — ) 

= E(ei^T^) 

= [ei^T^]-vK[(AT^TZ(AT^] 
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= [eisT(^-^]-\|/[sT(ALAT)s]. 
Thus,  Z is  Ep.^A  ju,  A Z A7). 


We  are  now  prepared  to  state  the  following  theorem. 


Theorem  2.2.6  Assume  the  observations  Yj, ...,  Yn  are  i.i.d.  from  a p-variate  elliptically 

symmetric  distribution.  Then,  under  H0,  defined  in  (2.1.4),  Vn,  defined  in  (2.1.5), 

computed  on  observations  Zj, ...,  Zn,  has  a small-sample  distribution-free  property  and  a 
2 

limiting  Xp^  distribution. 

Proof  of  Theorem  2.2.6  See  Randles  (1989),  p.  1046-1047. 


2,3  Asymptotic  Distribution  of  Vn  under  Contiguous  Alternatives 

In  this  section  we  will  find  the  asymptotic  distribution  of  Vn  under  a sequence  of 
alternatives  approaching  the  null  hypothesis  Hq  : 0 = 0.  In  doing  this,  we  will  restrict  our 
attention  to  a specific  class  of  elliptically  symmetric  distributions.  Let’s  assume  Y|, ...,  Yn 
are  i.i.d.  as  Y = (Yj, ...,  Yp)  , where  Y is  elliptically  symmetric  with  a density  function  fy 
of  the  form 

fy(y)  = Kplir1/2exp{-[(y  - u)7!'^  - u)/C0]v},  y e Rp,  (2.3.1) 


where 


pf(p/2v)  K vf(p/2) 

0 f[(p+2)/2v]  ’ p r(p/2v)(7tC0)p/2  ’ 


(2.3.2) 


r(w)  = j xw'1e'x  dx  for  w > 0, 
0 
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and  Rp  is  the  Euclidean  p-space.  It  can  be  verified  that  the  expression  in  (2.3.1)  is  a valid 
density  function  and  that  represents  the  mean  and  E,  the  variance-covariance  matrix.  This 
family  includes  the  multivariate  normal  distribution  (v  = 1),  heavier-tailed  distributions 
(0  < v < 1)  and  lighter-tailed  distributions  (v  > 1).  As  explained  in  the  previous  sections, 

we  will  need  to  derive  the  density  function  of  the  transformed  sample,  whose  form  will  be 
used  when  deriving  the  Pitman  asymptotic  relative  efficiency  of  Vn  relative  to  Hotelling’s 

T in  next  section. 


Lemma  2.3.3  For  the  family  of  distributions  given  in  expression  (2.3.1),  the  transformed 
sample  Zj, . . .,  Zn  has  a density  function  of  the  form 

fz(z)  = KpIA  I ATr1/2g[(z-Aii)T(A  Z A^fe-Au)],  z e Rp_1,  (2.3.4) 

where 

g(t)  = J exp{-[(t+s2)/C0]v}ds.  (2.3.5) 

-oo 


Proof  of  Lemma  2,3.3  Letting  x = 


then  we  can  write 


yi-yP 

yp-i-yp 

. yP  . 


y = 


.£l 


y = By, 


where  B = 


’A' 

.<£] 


T 

• £p  = (0,  ...,  0,  1)  , a pxl  vector  with  1 on  the  pth  component  and  0’s 


elsewhere,  and  A is  defined  in  (2.1.3).  Since  B is  nonsingular,  it  follows  that  y = B^x. 
The  jacobian  of  the  transformation,  denoted  by  J(y  =>  x),  is  equal  to  IB'1!  = IBP1  =1.  Hence 
the  density  function  of  X is 
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fx(x)  = fy(x)l  x=B_1x 

= KpISI'  1/2exp  { - [ (B ' 1 x-ja)TL'  1 (B ' 1 x-ji)/ C0] v } , 

= Kpl  Zl’  ^exp  { - [ (x- B yJT(B  E BT)' 1 (x-B  jjl)/ C0]  v } . 

Since  IB  E BTI  = IZI  due  to  IBI=  1,  it  follows  that 

fxOO  = KplB  Z BTI ' 1/2exp  { - [(x-B  ft)T(B  Z BT)‘ 1 (x-B  Ji)/Co] v } . 

T 

Letting  B Z B = V and  B = n,  we  rewrite  the  above  expression  as 

fx(x)  = KpIVr1/2exp{-[(x-Ti)TV-1(x-ii)/Co]v}. 

Thus,  the  density  function  of  Z is 
fz(2)  = | fx00  dyp 

-oo 

oo 

= Kpivr1/2  jexp{-[(x-n)TV-1(x-n)/Co]v}dyp.  (2.3.6) 

-OO 


Note  that 


V = BZBt  = 


AZAt  A Z ep 
epTIAT  eJZep 


and  n = Bu  = 


A ji 

L^p 


a 

.Mp. 


Denoting  V by 


Yu  Y 12 

Y21  V22 


, we  have  Vu  = A Z AT  is  nonsingular  and  V22  = e 1 Z e p 


-p- 


is  a positive  scalar.  Next  we  use  the  fact  that 


V1  = 


Yu  Y12 

Y21  V22  J 


-rl 


Y21 


22 


v;1!  + rl11Y12v-2‘2.1y21vj11  -Yi'vnV^., 


-viiY21rA 


y1 

v22-l 


where 

V22I  = V22~  Y2iY‘i\Yi2- 


We  can  expand  the  quadratic  form  in  expression  (2.3.6)  as 

(x-n)TY'1te-n) 

= (z-fl)TYA(s-fi)  + (l-fl)TY'1l1Yi2V2i!.1Y2iY'111(i-fi) 


+ (yp-h>)2v22-i  - ^-fflVAYiaV^yp-ty 

- + [(yP-Mp)v»r(rf>Tiiiii2^]2-  (2.3.7) 


Using  the  fact  that 

00  = Yi  ii-iv22  - y21v'11y12i  = iyuHV22.,i 


and  expression  (2.3.7),  we  may  write  expression  (2.3.6)  as 

fz(z) 

- KpIYul'1/2IV22.1l‘1/2 

• / expt-Ite-fflVAfe-fi)  + ((yp-Hp)vA'2,  - (2-ffiTY'111Yi2V212A)2)/C0]v)dyp. 

-oo 

Taking  s = (yp'M-p)V22^  - (z-fi)TV  /1Y12V224  > we  can  write  the  above  expression  as 
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oo 

fz(z)  = KplVnrw  j exp{-[((z-a)TV'1‘1(z-fi)  + s2)/C0]v)ds.  (2.3.8) 

-OO 

Since  Vj]  = ASA  and  0 = A ji,  expression  (2.3.8)  is  equivalent  to  expression  (2.3.4). 
This  completes  the  proof.  It  can  be  verified  that  expression  (2.3.8)  is  valid  density 
function. 

Now,  we  are  in  the  position  to  discuss  the  asymptotic  distribution  of  Vn  under  a 
sequence  of  alternatives  approaching  the  null  distribution.  Under  Hq,  Z\, ...,  Z„  are  i.i.d. 
with  density  of  the  form 

%(z)  = Kp  IAIATr1/2  ~ exp  { - [ (zt(A  Z At)'  1 z+s2)/C0] v } ds.  (2.3.9) 

-OO 

Under  a sequence  of  alternatives  let  Zj, ...,  be  i.i.d.  with  density  fztz-cn'1^),  where  fg 
is  given  in  (2.3.9)  and  c € Rp‘1-{0}  is  arbitrary,  but  fixed.  It  is  shown,  with  the  outline  of 
the  proof  given  in  Appendix  A,  that  if  4v  + p > 3,  then 

r cTdfz(z)  -,9 

0 < I£(fz)  = f L ~fz(z)  J fz(^)d2  < 00  for  all  c ^ 0. 

Thus,  when  4v  + p > 3,  the  rationale  of  Hajek  and  Sidak  (1967),  p.  212-213  and  earlier, 
shows  that  the  alternatives  are  contiguous  to  the  null  hypothesis.  Noting  that  both  Vn  and 

j 

Hotelling’s  T are  affine-invariant,  we  can,  without  loss  of  generality,  assume  throughout 
that  A ZAt  = Ip_i,  the  (p-l)x(p-l)  identity  matrix. 


Theorem  2.3. IQ  Let  Yj,  ...,  Yn  be  i.i.d.  with  density  function  fy  given  in  expression 
(2.3.1).  Assume  the  density  function  of  Zj,  ...,  Zn  has  the  form  of 

fz(z)  = Kpg(zTz)  = Kp  j exp  { -[(zTz+s2)/C0]V } ds,  z e Rp_1. 

-OO 
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If  4v  + p > 3,  then  under  the  sequence  of  contiguous  alternatives. 


d ? 

vn  — > Xp_i 


Rg'(R2) 

g(R2) 


(2.3.11) 


where 

g(R2)  = j exp{-[(R2+s2)/C0]v}ds, 

-oo 


and 


g'(R2)  = dg(R2)/dR2  and  R2  = ZTZ. 

Proof  of  Theorem  2.3.10  See  Randles  (1989),  p.  1050. 

Note  that  the  noncentrality  parameter  in  (2.3.1 1)  can  be  simplified  and  this  will  be 
done  in  the  next  section. 


2,4  The  Pitman  Asymptotic  Relative  Efficiency  of  Vn  Relative  to  Hotelling’s  T2 
In  this  section  we  will  use  Pitman  relative  efficiencies  to  make  comparisons 

9. 

between  Vn  and  Hotelling’s  T in  the  repeated  measurement  design  settings.  Because  these 

statistics  are  all  affine-invariant  we  may,  without  loss  of  generality,  make  the  simplifying 

assumption  that  the  transformed  sample  variance-covariance  matrix  is  the  identity.  We  will 

apply  the  asymptotic  results  under  the  contiguous  alternatives  in  Theorem  2.3. 10. 

The  Pitman  approach  to  asymptotic  relative  efficiency  compares  two  test  sequences 
(Snj)  and  {Tm.}  as  the  sequence  of  alternatives  H, : 0 = 0it  approaches  the  null,  which  we 

are  taken  to  be  H0  : 0 = 0.  The  subscripts  nj  and  mj  are  the  sample  sizes  for  tests  Sn.  and 

Tmj,  respectively.  Let  (3{Sn.,  0j } and  p{Tm.,  0 j } denote  the  powers  of  the  tests  based  on 
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Sn.  and  Tm.,  respectively,  when  0 = 0j.  Assume  that  nj  and  mj  are  such  that  the  two 
sequences  of  tests  have  the  same  limiting  significance  level  a and 

ac.lim  (3{Sn  , 0j}  = .lim  (3{Tm  , 0j]<l. 

1-400  1 1— >00  1 


Then  the  Pitman  asymptotic  relative  efficiency  (ARE)  of  {Sn.}  relative  to  {Tm.}  (or  simply 
of  S relative  to  T)  is 

ARE(S,  T)  = .lim  , 

1— >00  n j 

provided  the  limit  exists  and  is  the  same  for  all  such  sequences  {nj}  and  {nij},  and 
independent  of  the  {0,}  sequence.  (See,  e.g.,  Randles  and  Wolfe  [1979],  p.  144.) 

Hannan  (1956)  shows  that  if,  under  the  sequence  of  alternatives  Hj,  the  test 

sequences  { Sn. } and  [Tm.]  are  asymptotically  noncentral  Chi-square  with  the  same 

2 2 

degrees  of  freedom  and  noncentrality  parameters,  8S  and  5T,  respectively,  then 


ARE(S,  T)  = 


It’s  well  known  that  under  the  sequence  of  contiguous  alternatives  described  in  the  last 
section,  and  taking  AIAT  = Ip-i,  the  asymptotic  distribution  of  Hotelling’s  T2,  where 

T2  = n ZT  (£)_1Z  with  Z = n'1  £ Zj  and  £=  (n-1)’1  £ (Zj-Z)(Zj-Z)T,  is  noncentral 

i=l  i=l 

chi-square  with  p-1  degrees  of  freedom  and  noncentrality  parameter 


(2.4.1) 


(See,  e.g.,  Puri  and  Sen  [1971],  p.  173.)  To  derive  ARE(Vn,  T2)  we  use  the  following 


lemma. 
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Lemma  2.4.2  Let  Yj,  Yn  be  i.i.d.  with  density  function  given  in  expression  (2.3.1), 
and  R2  = ZTZ.  Taking  AIAT  = Ip_i,  then,  under  Hq  : 0 = 0,  the  density  function  of  R2  is 
of  the  form 

K tt(P'1)/2 

fR2(r)  =— f r^D^gCr),  r > 0,  (2.4.3) 

K r[(p-l)/2] 


where 

oo 

g(r)  = j exp  { -[(r+s2)/Q)] v } ds,  (2.4.4) 

-OO 

and  Kp  and  Co  are  defined  in  (2.3.2). 

Proof  of  Lemma  2.4.2  Taking  A E AT  = Ip-i  and  under  H0  : 0 = 0,  we  have  Z is 
Ep.i(Q,  Ip_i)  with  density  function  of  the  form 

%(z)  = Kpg(zTz)  = Kp  j exp  { - [ (zTz+s2)/CQ] v } ds,  z e Rp_1. 

-OO 

Thus,  R2  = ZTZ  has  density  of  the  form  given  in  expression  (2.4.3).  (See,  e.g.,  Muirhead 
[1982],  p.  37.) 


Theorem  2.4,5  Assume  Yj,  ...,  Yn  are  i.i.d.  from  a density  given  by  (2.3.1).  In  the 
repeated  measurement  settings,  taking  A E AT  = Ip.!,  if  4v  + p > 3,  then  the  Pitman 
asymptotic  relative  efficiency  of  Vn  relative  to  Hotelling’s  T2  is 


ARE(Vn,  T2)  = 


4r4(p/2)r2[(p-  i)/2v]r[(P+2)/2v] 
p(p-i)r4[(p-i)/2]r3(P/2v) 


(2.4.6) 


Proof  of  Theorem  2.4.5  If  4v  + p > 3,  taking  A E A.T  = Ip.!  and  using  expression 
(2.3.11),  we  have,  under  the  sequence  of  contiguous  alternatives,  the  test  Vn  is 
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asymptotically  noncentral  chi-square  with  p-1  degrees  of  freedom  and  noncentrality 
parameter  Sy  , where 


_±_ 

p-l 


Rg'(R2) 

g(R2) 


It  follows  that 


4, 


ARFfV  T^'l  - - 4 If  rRg'(R)l) 

ARE(Vn’T)‘  52  2 ■p-llVg(R2) 


(2.4.7) 


Under  Hq,  we  Find,  using  expressions  (2.4.3)  and  (2.4.4), 


= E 


r R-gW)  i 
EL  g(R2)  -I 

VrV(r2) 

g(R2) 


[ 


] 


°°  V r-g'(r) 

J— gW- fR2()dr 
VT-g'(r)  K-jifP-')/2 


i SCO 

K ttCp-D/2 


£2 r(P'1)/2‘1g(r)dr 

r((p-l)/2) 


r«p-l)/2)  )frp/2'1  g'(r)dr. 

0 


It  can  be  verified  that 

g'(r)  = ( J exp{-[(r+s2)/C0]v}ds)' 

-oo 


!°  r <KexP  { -[(r+s2)/C0]v } ) , 

= i[ s ]ds- 
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(See,  e.g.,  Trench  [1978],  Theorem  5.2  and  5.6  on  p.  581  and  p.  586,  respectively.) 
Letting 


1 E_ 

p r[(p-i)/2]  ’ 


(2.4.8) 


we  have 


R-g'(R2)  1 
g(R2)  J 


= Cp  J rP^1  | [ 

0 -OO 


1 °°r  0(exp(-[(r+s2)/Co]v}) 


3r 


] ds  dr 


OO  _ OO 


= Cp  j rP/2’1  j exp{-[(r+s2)/C0]v}(-v)((r+S2)/C0)v'1  C'Jds  dr 

0 -OO 


(_7  V r'  \ oooo 

— VS  S JrP/21  exP{-[(r+s2)/C0]v}(r+s2)v'1  ds  dr. 
C0  0 0 


Taking  r = r,  we  can  write 
,2 


E ['  R'g  (2R  } J j tp'J  exp { - [(t2+s2)/C0] v } (t2+s2)v_  1 ds  dt.  (2.4.9) 

g(R  ) Cq  oo 


Letting  s = VCq  x1/vsin(0)  and  t = ^[CQ  x1/vcos (0),  we  find  that 

3 s/00  = ^[CQ  x1/vcos(0),  0s/3x  = v_1VCo  x^1'v)/vsin(0)j 


0t/00  = -^Cq  x1/vsin(0),  0t/9x  = v'!VC o x^1"v^/vcos(0), 


and 


s2+t2  = C0  x2/v. 
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The  jacobian  of  the  transformation  is 

v'1Cox^2'v^/Vcos2(0)  + v'1Cox(2‘v^vsin2(0)  = v'1C0x^2'v^/v. 

Thus,  we  can  write  (2.4.9)  as 
r R'g^R2) 

EL  g(R2)  -I 

= ( 7^)j  j (VCoxl/Vcos(0))p'1(Cox2/v)v'1exp(-x2)v'1Cox(2'v)/vd0dx 

C0  0 0 

= (-4CpC(0P'1)/2  )•(  j x(p+v'1)/vexp(-x2)dx)  ( ]/cosp*1(0)d0).  (2.4.10) 

0 o 


The  constant  term  in  the  product  above  is 
-4CpC^1)/2 


= -4 


K^P-^2 

r((P-i)/2) 


p(p-l)/2 

U0 


(using  expression  (2.4.8)) 


1 vHp/2) 

r(P/2v)(7tc0)p/2 


jcG P^_ 
r[(P-i)/2] 


r(p-D/2 


(using  expression  (2.3.2)) 


1 vr(p/2)  _1  _ 

r(p/2v)T[(p-i)/2]  V K V Co 


The  first  integral  in  expression  (2.4.10)  is 

~(P+v4)/vexp(_x2)dx 

0 

= ~ f y(p"1)/2vexp(-y)dy  (using  y = x2) 

2 0 
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— j y^p+2v'1^2v"1exp(-y)dy 
2 0 


Xr(Pl2Yd)  = J— r(Eil 

2 2v  2 2v 


_l_(2d)r(2d)  = Mlr(Pd 

2 2v  2v  4v  2v 


The  second  integral  in  (2.4. 10)  is 


k/2  , jt/2.  , ^ n-f-) 

f cosP'H^Od©  = f sinP'^ejdO  = - — 

(j  J 2r<-^-) 


(See,  e.g.,  Beyer  [1987],  p.  289.)  Thus,  the  expected  value  in  expression  (2.4.10)  is 
r R-g'(R2)  ] 

Eh^rJ 

= r.4 vr<P/2> 1 

r(p/2v)r((p-l)/2)  V n V C0  4v  2v  J L 2T(-e±L)  J 


-(p-1)  r2ef-)  r(-^-> 


(2.4.11) 


Therefore,  using  expressions  (2.4.7)  and  (2.4.11),  we  have 


ARE(Vn,  T2) 


- -[ 

n-1  L 


(P-D  r2ef)  r(-^-) 


P‘1  2r(^)re|1-)reE2i)VQ 


] 


- [ 


(P-i)r4ef)r2(-ti-)  -.  r n-^2-) 


TH-]  • [ —~T~  ] • (using  (2.3.2)) 


r2^2^)^^)  pre£) 


2m 


Noting  that 
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_P±^  _ rr^i 


rev-)  = ne^-)+i]  = e 


— )re^ 


-), 


we  have 


ARE(Vn,T2) 


(P-D  r4ffo  r2(-^->  n-ffi 
pr3(^-)r4e^)(-ti-)2 

4 r4(-f->  r2(-E^-) 
p(p-i) 


This  ARE  is  evaluted  in  Table  2. 1 for  selected  values  of  v and  p satisfying  the  condition 
4v  + p > 3.  Note  that  when  p > 3,  4v  + p > 3 for  all  v > 0,  and  when  p = 2,  4v  + p > 3 for 
all  v > .25.  When  the  underlying  population  is  multivariate  normal  or  close  to  normal 
(v  = 1.0,  .75,  respectively),  Hotelling’s  T2  performs  better,  yet  Vn  appears  to  be  quite 

competitive  and  the  efficiencies  increase  as  p increases.  For  light-tailed  distributions 
(v  = 2,  3,  4,  and  5)  the  sign  test  is  not  as  effective  as  T2,  but  the  ARE’s  increase  with  p. 
For  heavy-tailed  distributions  (v  = .50,  .25,  .20,  .15,  and  .10)  Vn  is  more  effective  than 

T . Here  the  efficiencies  decrease  with  p.  In  fact,  it  can  be  verified,  using  Stirling’s 
formula  for  approximation  of  m!  and  the  result 

lim  (1+A.x)1/,X  = e^  for  constant  X, 
x— 

that  for  fixed  v,  ARE  — > 1 as  p — » For  fixed  p,  it  can  be  shown  that  as  v — > 0,  the 

ARE  ->  °°  and  as  v->  oo,  the  ARE  -»  {4p2r4(p/2)}/{(p+2)(p-l)3r4((p-l)/2)},  the  latter  is 
evaluated  for  p = 2 to  12  and  is  displayed  in  Table  2.1  under  the  column  v = oo. 
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*Conjectured  ARE's  when  p = 2 and  v < .25,  using  the  form  of  the  efficiency  given  in  Theorem  2.4.5. 


CHAPTER  3 

A MULTIVARIATE  SIGNED-RANK  TEST  BASED  ON  INTERDIRECTIONS 
FOR  REPEATED-MEASURES  DESIGNS 


3.1  Definition  of  the  Test  Statistic 

In  this  section  we  describe  the  multivariate  signed-rank  test,  denoted  by  Wn, 
proposed  by  Peters  and  Randles  (1990)  for  the  one-sample  multivariate  location  problem. 
As  explained  in  section  2.1,  this  test  can  also  be  applied  to  repeated  measures  designs  for 
detecting  treatment  effects. 

Using  the  same  notations  as  in  chapter  2,  we  let  Yj,  ...,  Y„  be  i.i.d.  from  a 
p-variate  elliptically  symmetric  distribution  and  Zj,  ...,  Zn  be  the  transformed  sample 
defined  in  (2.1.3).  Recalling  the  result  of  Theorem  2.2.5,  we  have  that  Zt,  ...,  Zn  are 

i.i.d.  from  a (p-l)-variate  elliptically  symmetric  distribution.  Thus,  it  is  logical  to  measure 
the  distance  of  each  observation  Zj,  i = 1,  ...,  n,  from  the  origin  in  terms  of  elliptical 

contours  and  to  use  the  ranks  of  these  distances  along  with  the  observations’  directions  in 
forming  a test  statistic. 

We  now  describe  such  a signed-rank  statistic  based  on  Vn  which  includes  the 

univariate  signed-rank  statistic  as  a special  case.  Specifically,  let  us  form  estimated 
Mahalanobis  distances  via 

Di  = ZT£'1Zi  ,i=l n,  (3.1.1) 

where 
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1 n j 


is  a consistent  estimator  of  the  null  hypothesis  variance-covariance  matrix  of  Zj,  provided  it 
exists,  with  H0  defined  in  (2.1.4).  Let  Qj  = Rank  (Dj)  among  £ji>  ....  Dn,  i = 1,  n. 
We  now  weight  the  (i,  k)th  term  in  the  sum  Vn,  defined  in  (2.1.5),  by  QjQk,  and  consider 
the  statistic 


n 


i=lk=l 


Icos(7tp.,  ) — 


Qi  Qk 


iky  n n 


(3.1.2) 


where  Pik  is  defined  in  (2.1.6).  We  reject  Hq  in  favor  of  Ha  for  large  values  of  the  statistic 
Wn. 

/\ 

Since  the  Pik,  i,  k = 1,  ...,  n,  are  invariant  with  respect  to  a nonsingular  linear 
transformation  (as  shown  by  Randles  [1989])  as  are  the  Di,  i = 1,  ....  n,  it  is  clear  that  Wn 
is  likewise  affine-invariant.  When  p = 2,  the  test  based  on  Wn  is  the  two-sided  univariate 
Wilcoxon  signed-rank  test.  For  p > 2,  Wn  does  not  have  a small-sample  distribution-free 

property,  but  its  large-sample  null  distribution  is  convenient,  as  is  shown  in  the  next 
section. 


3.2  Asymptotics 

In  this  section  we  develop  the  asymptotic  distributional  properties  of  nWn  under  the 

class  of  elliptically  symmetric  distributions.  Assume  that  Y is  elliptically  symmetric,  that  is, 
the  density  function  of  Y is  of  the  form 

fY(y)  = Kplir1/2h[(y-kL)TL'1(y-U)],  ye  Rp,  (3.2.1) 

where  j±  is  the  point  of  symmetry,  I is  the  variance-covariance  matrix  of  Y,  provided  it 
exists,  and  Kp  > 0.  The  null  distribution  of  nWn  is  established  in  the  following  lemma. 
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Theorem  3.2.2  Assume  the  observations  Yj, Yn  are  i.i.d.  from  a p-variate  elliptically 

symmetric  distribution  with  a density  function  defined  in  (3.2.1).  Then,  in  the  repeated 
measurement  settings,  the  test  nWn,  defined  in  (3.1.2),  has  a limiting  distribution 

under  Hq,  defined  in  (2.1.4). 

Proof  of  Theorem  3.2.2  Note  that,  under  H0,  taking  A Z AT  = lp,h  the  (p-l)x(p-l) 

identity  matrix,  Z\,  ...,  Z^  are  i.i.d.  as  Z = (Zj,  ...,  Zp.^),  where  Z is  from  an  elliptically 

symmetric  distribution  and  can  be  expressed  as  Z = RU,  where  R2,  as  before,  equals  ZTZ 

and  U is  distributed  uniformly  on  the  (p-l)-dimensional  unit-sphere  independent  of  R. 
(See,  e.g.,  Johnson  [1987].)  It  can  be  verified  that  E[(ZjZk)2]  = E[R4U2U2]  = 

E[R4]E[U2U2]  < E[R4]  < oo,  for  all  j,  k = 1,  ...,  p-1.  Thus,  via  the  Lindeberg-Levy 
Central  Limit  Theorem  (see,  e.g.,  Serfling  [1980],  p.  28),  each  element  of  Vn(£  - Ip^), 

^ 1 n -J>  /s 

where  £ = - X Z\  Zx , is  asymptotic  normal  under  Hq.  Therefore,  we  have  Vn(z  - Ip.i) 
= Op(l)  under  Hq.  (See,  e.g.,  Serfling  [1980],  p.  8.)  So,  the  test  nWn  has  a limiting  y2  . 

P-1 

distribution  under  Hq.  (See  Peters  and  Randles  [1990],  p.  553.) 

Next,  we  derive  the  Pitman  asymptotic  relative  efficiency  of  nWn  relative  to 
Hotelling’s  T . To  do  this,  we  first  establish  the  asymptotic  distribution  of  nWn  under 

contiguous  alternatives,  as  described  in  section  2.3,  for  a general  class  of  elliptically 
symmetric  distributions. 

Lgmma  3.2.3  Suppose  Xj,  ...,  Xn  are  i.i.d.  from  an  elliptically  symmetric  distribution 
with  density  of  the  form 


fx(x)  = Kmh[(x  - U)T(x  - u)],  x e Rm, 
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satisfying 


0 < I#x)  = 


dT2fx(x)  -|2 
%(*)  J 


fxOOdx  < oo 


for  all  de  Rm-{0}. 


(3.2.4) 


Under  the  sequence  of  contiguous  alternatives  for  which  X has  density  of  the  form 
fxfe-d/Vn),  we  have 


d 2 
nWn  ->  x2 


m 


( 12  fr,  r K(R2)Rh'(R2) 
VmlW  h(R2) 


]}Vd), 


(3.2.5) 


where 

H0  implies  U = Q,  R2  = XTX  , h'(R2)  = dh(R2)/dR2,  and  K(t)  = PHq(R2  < t),  t > 0. 


Proof  of  Lemma  3.2.3  Let 


i ; r dT2fx(Xi) 

L fx(Xi) 


]. 


* 


'/3S  [ 7=  £ Uj,K(R?) -f  £ UimK(R?)  ]T 

Vni=i  Vn  i=l  1 

~r=  X V3m  UjK(R?) , 

Vni=i  1 


and 


S„  - 


lT£, 


where  d,  ^ e Rm-{0},  and  U4  = [Uji,  ...,  Ujm]T  is  distributed  uniformly  on  the 
m-dimensional  unit-sphere  independent  of  Rv  Here  Xj  = RjU;  and  R2  = xJXj,  and  we 


can  write 
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n dTKmh'(x[xi)-2xi 


Tn  = tI 


[ 


^i=l  Kmh(x[xi) 


2h'(R?)RjdTU; 

= — z[ 1 

i=l  L h(R?) 


Thus,  under  Hq, 


T 

L 1 nJ 


= 7=2 

Vn  i=l 


V3m  K(R?)2LTUi 

2h'(Rf)RidTUi/h(R2) 


= ^|[wj  (“y). 


where 


Vi  = V3m  K(R?)^TUj  and  Wj  = 2h'(R|>RidTUi  /h(R?) 


2\  D jTi 


rVi 

Note  that,  under  Hq,  ^w. 


’s  are  i.i.d.  with 


:[wj  = 2 


and  Var 


Vj 

LWy 


an  a12 

a12  g22J 


where 


an  = E(Yf)  = 3mE[K2(Rf)]-EaTUiu[2J=2t12t 


since 


E[K2(Rf)]  = \ and  E(UiU?')  = i Imi 


022  = E(W2)=A{e[ 


h'(R?)-R, 


m L L h(R?) 


]} 


d*d, 


38 


and 


K(Rf)Rih'(Rf) 


h(R  f) 


Thus, 


all  a12 

L a12  a22- 


) under  Hq. 


(See,  e.g.,  Serfling  [1980],  Theorem  B,  p.  28.)  Applying  LeCam’s  Theorems  on 
contiguity  (see  Hajek  and  Sidak  [1967],  p.  212-213),  under  the  condition  defined  in 
(3.2.4),  we  have 


under  the  sequence  of  alternatives.  Therefore,  under  the  sequence  of  alternatives. 


(See  Peters  and  Randles  [1990],  Result  1 and  Theorem  2,  on  p.  555  and  p.  556, 
respectively.) 

Theorem  3.2.6  Assume  Yj,  ...,  Yn  are  i.i.d.  from  a density  given  by  (2.3.1).  In  the 
repeated  measurement  settings,  taking  A Z AT  = Ip.!,  if  4v  + p > 3,  then  the  Pitman 
asymptotic  relative  efficiency  of  nWn  relative  to  Hotelling’s  T2  is 


and  hence 
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16v  C2 

ARE(„Wn,  T2)  = ^ ( — - ) • 

co 

{OOOOOO  r ~ 

j j j f sP-^P-1  (r2+t2)v_1  exp{-[(s2+u2)/C0]v}exp{-[(r2 +t2)/C0]v}ds  du  dt  dr  } , 

0 0 0 0 

(3.2.7) 


where  Cp  is  defined  in  (2.4.8)  and  C0,  Kp  are  as  defined  in  (2.3.2). 

Proof  of  Theorem  3.2,6  Note  that,  under  H0,  taking  A I AT  = L.J,  Zh  ....  Zn  are  i.i.d. 
with  density  of  the  form 


oo 

f2(z)  = Kpg(zTz)  = Kp  j exp { -[(zTz+s2)/C0]v } ds,  z e RP'1.  (3.2.8) 

-OO 


(See  Lemma  2.3.3.)  Under  the  sequence  of  alternatives,  described  in  section  2.3, 
Zb  •••>  Zn  are  i.i.d.  with  density  %(z  - c/Vn  ),  where  c e RP'1-^}.  We  have  shown  that 

0 < Ic(%)  < 00  if  4v  + p > 3 (see  Appendix  A),  thus,  under  the  sequence  of  alternatives,  it 


follows  from  Lemma  3.2.3  that 


d 2 

nWn  “> 


(nfv  r K(R2)Rg'(R2) 
S-i1h0L  g(R2) 


where 


g'(R2)  = dg(R2)/dR2  and  R2  = ZTZ. 


Therefore, 

ARE(nWn,  T2)  = { E„  [—  ^r2)R§/(r2)  1 } 2 

P-1  H0  g(R2)  JJ 

Recall  Lemma  2.4.2  that  R2  = ZTZ  has  density  of  the  form 


(3.2.9) 
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K te^P-1)/2 

fR2(r)  =— j* r > 0. 

K F[(p-l)/2] 


Thus,  we  have 


„ r K(R2)Rg'(R2)  1 
W g(R2)  J 


J°  [ K(r)Vr  g'(r)  j 


fR2(r)  dr 


r(p-l)/2-lg(r)  ^ 


= Cp  j°°  KWg'WrP/2-1  dr. 
0 


Note  that 


g'(r)  = ( j (r+t2)v_1  exp{-[(r+t2)/Co]v}  dt. 

co  0 


It  follows  that 


[ K(R2)Rg/(R2)  j 


(-2v-C  \ 0000 

j j K(r)-rP/2-1(r+t2)v-1exp{-[(r+t2)/C0]V}  dt  dr 

C0  0 0 

(_2v  • C \ 0000  r 

J / [j  fR2(s)  ds  ] rp/2'1(r+t2)v'1  exp{-[(r+t2)/C0] 
<-0  0 0 0 


(3.2.10) 


} dt  dr. 


Using  the  expressions  (3.2.10)  and  (3.2.8)  for  fR2  and  g,  respectively,  in  the  above 
expression  yields 


41 


E r K(R2)Rg/(R2)  1 

H°L  g(R2)  J 


/ -4vCp\  -oooor  , 

= ( — ) / J jf  s(p-1)/2A  rp/2'1  (r+t2)v'1  exp  { - [(s+u2)/C0] v } • 

C0  0 0 0 0 


exp{-[(r+t2)/Co]v}ds  du  dt  dr 

, ■>6vCL 

= ( — ) j j j j sp-2rp-‘  (r2 +t2)v"'  exp{-((s2+u2)/C0]v)- 

c0  oooo 


exp{-[(r2+t2)/Co]V}  ds  du  dt  dr. 


Substituting  the  above  expression  into  (3.2.9)  yields  Theorem  3.2.6. 


3.3  Numerical  Evaluation  of  ARE(nWn.  T2) 

In  this  section  we  describe  the  numerical  evaluation  of  the  ARE  expression  given  in 
Theorem  3.2.6.  All  calculations  were  performed  with  an  IBM  computer  running  on  a 
VM/CMS  operating  system  using  fortran  77.  Simpson’s  rule  and  the  IMSL  subroutine 
DMLIN  were  used  to  integrate  single  integrals  and  3-dimensional  integrals,  respectively. 

Now  we  describe  how  the  calculations  were  performed  as  well  as  the  error 
associated  with  them.  Letting 

* -U  * -t  i * -r 

u = e , t = e , and  r = e , 

we  can  write 


VARE(nWn,T2) 
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/ 12  1 

r p 

/ P-1  ' 

1 cv 
^0 

){  ///  j”p-2(-Cnr*)p'1(r*t*u*)-,[(/rt’)2+(Cnr*)2iv-‘- 
0 0 0 0 


exp{-[(s2+(rmi*)2)/C0]v}exp{-[((rnr*)2+(^it*)2)/Co]v}ds  du*dt  *dr*  } 

s j j f [ | 4Cpsp'V1exp{-[(s2+(rnu)2)/C0]v}  ds  ]• 

0 0 0 0 

[V^^(-^)P'1[(^)2+(^)2]V‘1(t  r)'1exp{-[((M2+(/'nt)2)/C0]v}]du  dt  dr 
C0 


1 1 1 -(nr 

j j J [ j k(u,  s)  dsj  h(t,  r)  du  dt  dr, 
0 0 0 0 


(3.3.1) 


where 


k(u,  s)  = 4Cpsp‘2u'1exp{-[(s2+(/nu)2)/C0]v}  , 


h(t’ r)  = VS  (-/«r)p'1[(/nt)2+(/nr)2]v’1(t  r)-1exp{-[((r«r)2+(rnt)2)/C0]v}, 


4vC 

p-i  rv 


and  (n  is  the  natural  logarithm.  Note  that,  following  standard  arguments  and  using 
expressions  (2.4.7)  and  (2.4.9), 

1 1 

ss 

0 0 


/ /h(t,  r)  dt  dr  = V 3ARE(Vn,  T2) , 


(3.3.2) 


where  ARE(Vn,  T ) is  given  in  Theorem  2.4.5.  Using  Simpson’s  rule  with  an  absolute 
error  3ARE(Vn,  T2),  we  can  approximate  the  inner  single  integral  in  expression 
(3.3.1)  by  g(u,  r)  satisfying 


-(nr 


J 


k(u,  s)  ds  - g(u,  r)  I < ej/^J  3ARE(Vn,  T2). 


(3.3.3) 


43 


Hence,  using  expressions  (3.3.2)  and  (3.3.3),  we  can  write  expression  (3.3.1)  as 
VARE(nWn,  T2) 

- J j f [ g(u>  r)  ± ei/V  3ARE(Vn,  T2)  ] h(t,  r)  du  dt  dr 
0 0 0 

= jji  g(u,  r)h(t,  r)  du  dt  dr  ± [£l/V  3ARE(Vn,  T2)  ] j 1 /h(t,  r)  dt  dr 
0 0 0 0 0 

1 1 1 

= j j j g(u,  r)h(t,  r)  du  dt  dr  ± eh  (3.3.4) 

0 0 0 

Now,  using  the  IMSL  subroutine  DMLIN  with  an  absolute  error  £2,  we  can  approximate 
the  above  3-dimensional  integral  by  VAL  satisfying 

I 111  1 

I J j | g(u,  r)-h(t,  r)  du  dt  dr  - VAL  I < e2.  (3.3.5) 

0 0 0 

Combining  expressions  (3.3.4)  and  (3.3.5),  we  have 

V ARE(nWn,  T2)  = VAL  ± (e^), 

where  VAL  satisfies  expression  (3.3.5).  Thus,  we  approximate  ARE(nWn,  T2)  by  (VAL)2 
with  the  maximum  error 

ERREST  = 2VAL(E]+£2)  "h  (£i+£2)^-  (3.3.6) 

In  Table  3.1  we  present  the  asymptotic  relative  efficiencies  for  nWn  with  respect  to 

2 

T for  selected  values  of  v and  p satisfying  the  condition  4v+p  > 3.  We  also  include 
ARE(Vn,  T ),  in  parentheses,  for  easy  comparison.  When  the  underlying  population  is 
multivariate  normal  (v  = 1.0),  Hotelling’s  T2  performs  well.  Both  the  sign  test  and 
signed-rank  test  appear  to  be  quite  competitive,  with  the  signed-rank  test  slightly  better  than 
the  sign  test  provided  the  dimension  is  not  too  large.  For  light-tailed  distributions 
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(v  = 2.0  and  5.0)  and  p > 2,  the  signed-rank  test  Wn  has  greater  power  compared  to  both 
Hotelling’s  T and  Vn.  For  heavy-tailed  distributions  (v  = .50  and  .10),  Vn  is  clearly  most 
powerful,  although,  Wn  still  performs  well  relative  to  Hotelling’s  T2,  provided  the 
dimension  is  not  too  large. 

In  Table  3.2  we  display  the  values  of  ERREST,  defined  in  expression  (3.3.6), 
which  bound  the  error  in  the  estimates  of  ARE(nWn,  T2).  Assuming  Ej  = £2  = £,  where 
El  and  £2  are  defined  in  (3.3.3)  and  (3.3.5),  respectively,  we  use  the  foiling  e values  : 
E = .001  when  v = 1.0,  £ = .001  when  v - 2.0  and  p =2  to  6,  £ = .005  when  v = 2.0  and 
p =7  or  8,  £ = .01  when  v = 5.0  and  p = 2 to  6,  £ = .02  when  v = 5.0  and  p = 7 or  8, 
£ = .001  when  v = .50,  £ = .10  when  v = .10  and  p = 2,  and  £ = .015  when  v = .10  and 
p = 3 to  8. 


Table  3.1 

ARE(nWn,  T^),  with  ARE(Vn,  T^)  in  parentheses. 
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Table  3.2 

Error  Estimate  ERREST  of  ARE(nWn,  T2) 
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CHAPTER  4 
MONTE  CARLO  STUDY 


In  this  chapter  we  display  results  from  a Monte  Carol  study  when  the  dimension  is 
p = 4,  the  sample  size  is  n = 20,  and  the  significance  level  is  a = .05.  In  addition  to  the 
affine-invariant  statistics  T , Vn,  and  nWn,  we  examined  two  other  nonaffine-invariant 

statistics.  The  ANOVA  F test  is  included  along  with  the  rank  transformation  test  RT,  which 
were  both  introduced  in  chapter  1.  These  five  test  statistics  were  compared  under  five 
different  distributions.  They  were  quadrivariate  normal  distribution,  elliptically  symmetric 
distributions  with  density  of  the  form  given  in  (2.3.1),  Pearson  Type  II  and  Type  VII 
(see  Johnson  [1987],  p.lll  and  p.  1 17-1 18,  respectively),  and  the  mixtures  of 
quadrivariate  normal  distributions.  (A  quadrivariate  normal  mixture  is  obtained  by  selecting 
randomly  one  of  two  quadrivariate  normal  distributions.  Each  of  the  observations  is 
sampled  with  probability  p from  the  first  distribution  and  with  probability  1-p  from  the 
second  distribution.  See  Johnson  [1987],  p.56-57.)  These  distributions  were  located  at 
Ifa  = (m8,  m5,  m8,  0),  m = 0,  1,  2,  3,  for  the  original  sample  Yj,  ...,  Yn,  on  which 
ANOVA  F test  and  rank  transformation  test  RT  were  applied,  and  they  were  located  at  63  = 
(m8,  m8,  m8),  m = 0,  1,  2,  3,  for  the  transformed  sample  Zlf  ...,  Zn,  on  which  the  tests 
T , Vn,  and  nWn  were  performed.  The  value  of  8 was  adjusted  for  different  distributions  to 
examine  somewhat  silimar  points  on  the  power  curves.  Since  the  performances  of  tests 
ANOVA  F and  RT  depend  on  the  variance-covariance  structure  of  the  distribution,  for  each 
of  the  above  first  four  distributions,  we  considered  two  types  of  variance-covariance 
matrices,  one  with  I4,  the  identity  matrix,  and  the  other  one  with  a non  H-type  structure. 

For  the  mixtures  of  normal  distributions,  we  consider  one  mixture  with  a non  H-type 
structure,  and  three  other  mixtures  with  H-type  structures.  In  each  Monte  Carlo  simulation, 
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the  proportion  of  times  out  of  1000  in  which  each  test  statistic  exceeded  the  upper 
a-percentile  of  its  null  distribution  is  reported.  The  asymptotic  null  distribution  is  used 


to  determine  the  critical  value  for  tests  Vn  and  nWn.  For  tests  ANOVA  F and  rank 

transformation  RT,  the  null  distribution  Fj^y  is  used  to  determine  the  critical  value.  While 
2 57 

for  Hotelling’s  T , we  use  the  null  distribution  — F3 ,17  (namely,  a multiple  of  F3,j7)  to 


determine  the  critical  value.  All  calculations  and  random  variables  generations  were 
performed  with  an  IBM  computer  running  on  a VM/CMS  operating  system  using  fortran 
77.  Several  IMSL  subroutines,  to  be  described  later  in  this  chapter,  were  used. 

In  Tables  4.1  and  4.2,  the  results  from  the  Monte  Carol  studies  for  quadrivariate 
normal  distributions  with  Z = I 4 and  Z = E ET,  respectively,  are  presented.  Note  that  we 

used 


E = 


r 1 3 0 

-2  1 4 

1-1  1 
L -2  4 -1 


2 1 

-1 

0 ’ 
3 J 


(4.1.1) 


producing 


r 14  -1 

-1  22 
-2  1 
. 16  1 


-2  16  -1 
1 1 
3 -7 

-7  30  J 


(4.1.2) 


which  is  non  H-type.  The  IMSL  subroutine  GGNML  was  used  to  generate  N(0,  1) 
variables.  As  indicated,  1000  samples  of  size  20  and  significance  level  .05  were  used. 
When  Z = I4,  as  we  may  expect,  ANOVA  F and  rank  transformation  test  RT  have  better 
power  than  Hotelling’s  T2.  And  they  all  perform  better  than  Vn  and  nWn  (see  Table  4.1). 
However,  when  Z = E ET,  Hotelling’s  T2  has  good  power,  followed  by  Vn  and  nWn,  both 
performing  better  than  F and  RT  (see  Table  4.2).  This  illustrates  the  strong  dependence  of 
the  performances  of  tests  F and  RT  on  the  variance-covariance  structure  of  the  distribution. 
For  both  variance-covariance  structures,  ANOVA  F performs  better  than  RT,  and  nWn 


49 


performs  slightly  better  than  Vn,  which  agrees  with  the  result  found  in  Table  3.1  for  v =1.0 
and  p = 4. 

In  Tables  4.3  and  4.4,  we  display  the  Monte  Carlo  results  for  five  members  of  the 
class  of  elliptically  symmetric  distributions  with  density  of  the  form  given  in  (2.3.1)  with 
Z=l4  andZ=EE  , respectively.  Heavy-tailed  distributions,  v = .10  and  .50,  the 

quadrivariate  normal  distribution,  v = 1.0,  (repeated  in  this  table  to  permit  comparisons), 
and  light-tailed  distributions,  v = 50.0  and  100.0,  were  included.  Taking  y = 0 and  Z = I4 
in  (2.3.1),  it  follows  from  standard  arguments  that 

P(R2  < w)  = P(C0G1/v  < w),  w > 0, 

where  R2  = YT  Y,  Cq  is  defined  in  (2.3.2),  and  G has  the  distribution  of  Gamma(l,  p/2v). 

Thus  Y can  be  generated  via  Y = R U , where  R is  independent  of  U having  the  distribution 
of  a/CqG1/2v,  and  U is  uniformly  distributed  on  the  unit-sphere  Sp  (see  Johnson  [1987], 

p.  1 10).  Here,  the  IMSL  subroutines  GGAMR  and  GGSPH  were  used  to  generate 
Gamma(l,  r)  and  the  variables  that  are  uniformly  distributed  on  the  unit- sphere  Sp, 

respectively.  Although  it  is  difficult  to  compare  powers  for  statistics  with  different  rejection 
proportions  under  the  null  hypothesis,  it  appears  that  the  Monte  Carol  results  tend  to  agree 
with  the  asymptotic  results  of  the  previous  two  chapters.  For  light-tailed  distributions  (v  = 
50.0  and  100.0),  ANOVA  F and  RT  have  great  power  if  Z = I4,  while  Hotelling’s  T2  and 
nWn  have  the  greatest  power  when  Z = EET.  For  the  heaviest-tailed  distribution  (v  = .10) 
examined,  Vn  performs  best,  followed  in  order  by  RT  and  Hotelling’s  T2  when  Z = I4,  and 
by  Hotelling’s  T2  and  RT  when  Z = EET.  Generally  speaking,  for  light-tailed  distributions 
(v  = 50.0  and  100.0),  the  signed-rank  test  nWn  performs  better  than  the  sign  test  Vn,  and 
ANOVA  F performs  better  than  RT;  however,  for  heavy-tailed  distributions  (v  = .10  and 
.50),  the  interdirection  sign  test  Vn  performs  better  than  nWn,  and  RT  works  better  than  F. 
The  superiority  of  tests  ANOVA  F and  nWn  (RT  and  Vn)  for  very  light-tailed  (heavy-tailed) 
distributions  is  shown  here.  Except  for  the  heaviest-tailed  distribution  (v  = .10),  ANOVA  F 


50 


and  RT  both  perform  better  than  Vn  and  nWn  when  Z = I4,  however,  Vn  and  nWn  both 
perform  better  than  ANOVA  F and  RT  when  Z = EE  . Also,  we  notice  that  RT  performs 
better  than  Hotelling’s  T2  only  when  Z = I4. 

Next,  we  consider  the  light-tailed  Pearson  Type  II  distributions,  a special  case  of 
elliptically  symmetric  distributions,  with  the  density  function 

fY(x)  = r(m+3)JZr1/2(i-(y-u)Tr1(y-n)}m,  (4.1.3) 

— r(m+l)jr 

having  support  ^ 1 and  shape  parameter  m > -1.  The  Monte  Carlo  results 

for  this  family  of  distributions  with  Z = I4  and  Z = E ET  are  presented  in  Tables  4.5  and 
4.6,  respectively.  Taking  \x  = 0 and  Z = 1^  in  (4.1.3),  it  can  be  shown  that  R2  = YT  Y has 
the  distribution  of  Beta(2,  m+1).  Thus,  the  variable  generation  is  similar  to  that  in  the  last 
case  considered.  (See,  Johnson  [1987],  p.  1 16.)  Here,  the  IMSL  subroutine  GGBTR  was 
used  to  generate  beta  variables.  An  examination  of  Tables  4.5  and  4.6  indicates  that 
ANOVA  F and  RT  have  great  power  if  Z = I4,  while  Hotelling’s  T2  and  nWn  have  the 
greatest  power  when  Z = E E . The  signed-rank  test  nWn  performs  better  than  the  sign  test 
Vn,  and  ANOVA  F performs  better  than  RT.  Furthermore,  ANOVA  F and  RT  both 
perform  better  than  Vn  and  nWn  when  Z = I4,  while  Vn  and  nWn  both  perform  better  than 
ANOVA  F and  RT  when  Z = E ET.  These  results  agree  with  the  findings  in  Tables  4.3  and 
4.4  for  the  light-tailed  elliptically  symmetric  distributions. 

For  the  heavy-tailed  Pearson  Type  VII  distributions,  the  Monte  Carlo  results  are 
presented  in  Tables  4.7  and  4.8  for  Z = I4  and  Z = EE^,  respectively.  The  density 

function  has  the  form  of 

fY(*)  = — ~ 2'a'1/2l  1+  (x-u)TS'1(y-u))'m,  (4.1.4) 

— r(m-2)7T 

where  the  shape  parameter  m > 2.  Note  that,  taking  ji  = 0 and  Z = I4  in  (4.1.4),  Y can  be 
generated  via  Y = Z/Vs" (see  Johnson  [1987],  p.  118),  where  Z is  quadrivariate  normal 
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with  mean  0 and  variance-covariance  matrix  I4,  and  S is  independent  of  Z having  the 
distribution  of  %b,  b = 2m  - 4.  Note  that  when  m = 2.5,  Y has  a multivariate  cauchy 

distribution.  In  Tables  4.7  and  4.8,  the  sign  test  Vn  has  the  greatest  power.  Notice  that  Vn 
performs  better  than  nWn,  and  RT  works  better  than  ANOVA  F.  Also,  Vn  and  nWn  both 
perform  better  than  F and  RT  when  Z = E ET.  Note  that  the  same  general  trends  as  in  the 
heavy-tailed  elliptically  symmetric  distributions  are  seen  here. 

Finally  in  Table  4.9  we  examine  samples  from  quadrivariate  normal  mixtures 
violating  the  assumption  of  elliptical  symmetry  required  for  the  asymptotic  results  of  the 
previous  two  chapters.  As  we  may  expect,  the  sign  test  Vn  and  signed-rank  test  nWn 

appear  to  be  quite  robust.  For  each  mixture  in  Table  4.9,  the  parameters  of  the  distributions 
(Ul»  U2>  £i»  a°d  la)  and  the  mixing  probability  (p)  as  well  as  the  parameters  of  the  mixture 
(ji_and  Z ),  computed  by 

H = PUl  + (1-P)U2> 
and 

Z = pZt  + (l-p)Z2  + p(l-p)(Ui  - U2)(Ui  - U2)T,  (4.1.5) 

(see,  Johnson  [1987],  p.57),  are  indicated.  In  mixture  1,  a light-tailed  distribution  with 
H-type  structure,  we  see  that  tests  ANOVA  F,  RT,  and  Hotelling’s  T2  have  great  power, 
and  all  performing  better  than  Vn  and  nWn.  This  agrees  with  that  of  Table  4.1,  for  normal 
distribution  with  Z = I4.  An  examination  of  the  table  for  mixtures  2 and  3,  heavy-tailed 

distributions  with  H-type  structure,  shows  the  same  general  trends.  That  is,  the  superiority 
of  tests  Vn  and  RT  over  the  signed-rank  test  nWn  and  the  other  competitors  are  seen.  Also, 
by  comparing  with  the  result  for  elliptically  symmetric  distribution  with  v = .10  and  Z = I4 
in  Table  4.3,  and  that  for  Pearson  Type  VII  distribution  with  m = 3.0  and  Z = I4  in  Table 
4.7,  it  illustrates  that  the  performances  of  Hotelling’s  T2  and  ANOVA  F are  affected  by  the 
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asymmetry  of  the  distributions,  which  as  a result  makes  nWn  a better  test  than  Hotelling’s 
2 

T and  ANOVA  F.  Finally,  in  mixture  4,  a heavy-tailed  distribution  with  non  H-type 
structure,  the  robustness  and  superiority  of  tests  Vn  and  nWn  is  easily  seen. 

In  summary,  for  normal  to  light-tailed  symmetric  distributions  and  a variance- 
covariance  structure  that  is  non  H-type,  Hotelling’s  T2  has  the  greatest  power,  followed  by 
the  signed-rank  test  nWn.  When  the  structure  is  H-type  and  the  populations  are  normal  or 

light-tailed  the  best  is  ANOVA  F,  followed  by  RT.  For  the  heavy-tailed  distributions  with 
H-type  variance-covariance  structure,  the  interdirection  sign  test  Vn  has  the  greatest  power, 

followed  by  RT.  For  heavy-tailed  symmetric  distributions  with  non  H-type  variance- 
covariance  structure,  the  interdirection  sign  test  Vn  still  has  the  greatest  power,  but  this  time 

Hotelling’s  T is  second  best.  For  the  asymmetric  distribution  with  non  H-type  variance- 
covariance  structure,  the  interdirection  sign  test  Vn  and  the  signed-rank  test  nWn  are  the 
best.  Generally,  the  tests  Vn  and  nWn  both  perform  better  than  ANOVA  F and  RT  when 

the  variance-covariance  structure  is  non  H-type  while  the  contrary  is  true  only  for  normal  to 
light-tailed  distributions  with  H-type  variance-covariance  structures.  This  suggests  that  Vn 
and  nWn  are  promising  procedures  for  the  repeated  measures  settings.  The  tests  ANOVA  F 
and  RT  follow  the  same  general  patterns  as  nWn  and  Vn,  respectively.  That  is,  the  RT  (Vn) 
test  performs  better  than  ANOVA  F (nWn)  for  heavy-tailed  distributions,  and  ANOVA  F 
(nWn)  performs  better  than  RT  (Vn)  for  light-tailed  distributions.  Also  note  that  test  RT 

works  better  than  Hotelling’s  T2  for  distributions  with  H-type  variance-covariance 
structures.  Finally,  we  should  point  out  that  the  performance  of  the  test  based  on  nWn  is 
disappointing,  since  the  power  of  nWn  is  not  always  better  than  that  of  Hotelling’s  T2  for 
light-tailed  symmetric  distributions  (see  Tables  4.3  and  4.4).  One  explaination  to  this  may 
be  because  of  the  higher  probability  of  Type  I error  of  Hotelling’s  T2  test. 
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Table  4.1 

Monte  Carlo  Results  for  Quadrivariate  Normal  Distribution  with  £ = I4. 


Amount 

of 

Shift 

Statistics 

F 

RT 

»p2 

vn 

nWn 

.00 

.049 

.046 

.056 

.052 

.058 

.25 

.106 

.113 

.108 

.082 

.108 

.50 

.335 

.321 

.300 

.219 

.247 

.75 

.649 

.644 

.579 

.457 

.485 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.2 

Monte  Carlo  Results  for  Quadrivariate  Normal  Distribution  with  Z = E ET. 


Amount 

of 

Shift 

Statistics 

F 

RT 

T2 

vn 

nWn 

.00 

.081 

.092 

.055 

.040 

.056 

.60 

.091 

.101 

.108 

.088 

.114 

1.20 

.164 

.144 

.283 

.204 

.235 

1.80 

.279 

.248 

.540 

.430 

.453 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 


i 
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Table  4.3 

Monte  Carlo  Results  for  Elliptically  Symmetric  Distributions  with  Density  of  the 
Form  Given  in  (2.3.1)  with  Z = I4. 


Amount 

of 

Shift 

Statistics 

F 

RT 

T2 

vn 

nWn 

v = .10 

.00 

.026 

.048 

.016 

.032 

.040 

.17 

.079 

.226 

.097 

.276 

.124 

.34 

.248 

.636 

.367 

.669 

.288 

.51 

.500 

.872 

.626 

.894 

.432 

< 

11 

Ln 

0 

.00 

.053 

.057 

.050 

.032 

.048 

.25 

.110 

.119 

.101 

.079 

.088 

.50 

.321 

.368 

.303 

.260 

.234 

.75 

.684 

.699 

.618 

.573 

.462 

v = 1.0 

.00 

.049 

.046 

.056 

.052 

.058 

.25 

.106 

.113 

.108 

.082 

.108 

.50 

.335 

.321 

.300 

.219 

.247 

.75 

.649 

.644 

.579 

.457 

.485 
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Table  4.3  - - continued. 


Amount 

of 

Shift 

Statistics 

F 

RT 

rp2 

vn 

nWn 

v = 50.0 

.00 

.046 

.049 

.059 

.032 

.044 

.25 

.100 

.085 

.092 

.057 

.087 

.50 

.316 

.279 

.247 

.166 

.261 

.75 

.656 

.594 

.559 

.361 

.568 

v = 100.0 

.00 

.050 

.053 

.053 

.032 

.048 

.25 

.095 

.089 

.091 

.056 

.088 

.50 

.317 

.282 

.256 

.159 

.262 

.75 

.650 

.584 

.573 

.354 

.557 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.4 

Monte  Carlo  Results  for  Elliptically  Symmetric  Distributions  with  Density  of  the 
Form  Given  in  (2.3.1)  with  £ = E ET. 


Amount 

of 

Shift 

Statistics 

F 

RT 

.p2 

vn 

nWn 

v = .10 

.00 

.055 

.081 

.021 

.049 

.037 

.40 

.073 

-.128 

.094 

.277 

.127 

.80 

.124 

.277 

.325 

.657 

.268 

1.20 

.224 

.462 

.586 

.879 

.418 

v = .50 

.00 

.079 

.081 

.060 

.049 

.063 

.60 

.094 

.104 

.110 

.100 

.107 

1.20 

.154 

.152 

.300 

.266 

.238 

1.80 

.268 

.219 

.609 

.558 

.452 

v = 1.0 

.00 

.081 

.092 

.055 

.040 

.056 

.60 

.091 

.101 

.108 

.088 

.114 

1.20 

.164 

.144 

.283 

.204 

.235 

1.80 

.279 

.248 

.540 

.430 

.453 
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Table  4.4  - - continued. 


Amount 

of 

Shift 

Statistics 

F 

RT 

T2 

vn 

nWn 

v = 50.0 

.00 

.078 

.079 

.060 

.049 

.048 

.80 

.110 

.105 

.139 

.098 

.124 

1.60 

.197 

.163 

.438 

.270 

.434 

2.40 

.390 

.275 

.833 

.620 

.830 

v = 100.0 

.00 

.073 

.080 

.065 

.049 

.055 

.80 

.106 

.109 

.138 

.099 

.123 

1.60 

.203 

.157 

.428 

.275 

.447 

2.40 

.394 

.280 

.824 

.617 

.822 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.5 

Monte  Carlo  Results  for  Pearson  Type  13  Distribution  with  Shape  Parameter  m,  and 
Density  of  the  Form  Given  in  (4.1.3)  with  £ = 14. 


Amount 

of 

Shift 

Statistics 

F 

RT 

rji2 

vn 

nWn 

m = 1.0 

.00 

.045 

.045 

.050 

.032 

.039 

.10 

.115 

.106 

.108 

.064 

.092 

.20 

.402 

.345 

.340 

.221 

.292 

.30 

.771 

.724 

.687 

.493 

.649 

m = 2.0 

.00 

.052 

.045 

.051 

.032 

.053 

.10 

.136 

.130 

.122 

.077 

.113 

.20 

.498 

.451 

.415 

.281 

.370 

.30 

.867 

.830 

.788 

.601 

.726 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.6 

Monte  Carlo  Results  for  Pearson  Type  II  Distribution  with  Shape  Parameter  m,  and 
Density  of  the  Form  Given  in  (4.1.3)  with  £ = E ET. 


Amount 

of 

Shift 

Statistics 

F 

RT 

rji2 

v„ 

nWn 

m = 1.0 

.00 

.078 

.080 

.060 

.049 

.053 

.30 

.111 

.109 

.153 

.111 

.144 

.60 

.223 

.170 

.490 

.325 

.454 

.90 

.438 

.330 

.877 

.704 

.832 

m = 2.0 

.00 

.072 

.075 

.059 

.049 

.054 

.30 

.118 

.113 

.172 

.125 

.160 

.60 

.275 

.195 

.587 

.425 

.522 

.90 

.540 

.403 

.930 

.812 

.885 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.7 

Monte  Carlo  Results  for  Pearson  Type  VII  Distribution  with  Shape  Parameter  m,  and 
Density  of  the  Form  Given  in  (4.1.4)  with  Z = I4. 


Amount 

of 

Shift 

Statistics 

F 

RT 

T2 

vn 

nWn 

m = 2.5 

.00 

.019 

.052 

.022 

.052 

.050 

.40 

.024 

.106 

.042 

.107 

.097 

.80 

.067 

.262 

.109 

.287 

.190 

1.20 

.120 

.476 

.214 

.520 

.312 

m = 3.0 

.00 

.035 

.046 

.030 

.052 

.058 

.50 

.179 

.324 

.226 

.295 

.220 

1.00 

.592 

.863 

.661 

.829 

.575 

1.50 

.851 

.985 

.901 

.976 

.777 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.8 

Monte  Carlo  Results  for  Pearson  Type  VII  Distribution  with  Shape  Parameter  m,  and 
Density  of  the  Form  Given  in  (4.1.4)  with  Z = EET. 


Amount 

of 

Shift 

Statistics 

F 

RT 

T2 

vn 

nWn 

m = 2.5 

.00 

.025 

.077 

.027 

.040 

.048 

1.20 

.031 

.098 

.058 

.135 

.110 

2.40 

.054 

.173 

.149 

.376 

.232 

3.60 

.087 

.294 

.281 

.653 

.375 

m = 3.0 

.00 

.059 

.080 

.027 

.040 

.051 

1.10 

.093 

.137 

.182 

.251 

.199 

2.20 

.228 

.340 

.567 

.709 

.508 

3.30 

.442 

.620 

.848 

.939 

.715 

Entries  : Proportion  of  times  each  test  statistic  exceeded  the  upper  a-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 

Significance  Level : a = 0.5. 
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Table  4.9 

Monte  Carlo  Results  for  Quadrivariate  Normal  Mixtures  with  Mean  y,  and  Variance- 
Covariance  Matrix  Z given  in  (4.1.5). 


Amount 

of 

Shift 

Statistics 

F 

RT 

»p2 

vn 

nW„ 

Mixture  1 

P = -5,  JA]  = [1,  1,  1,  1]T  = -y.2,  Zi  = £2  = 14. 

r 2 1 1 ll 

* --  112  1- 
.1112. 

.00 

.049 

.051 

.047 

.033 

.049 

.35 

.190 

.179 

.171 

.121 

.160 

.70 

.604 

.565 

.528 

.405 

.428 

1.05 

.929 

.915 

.896 

.800 

.770 

Mixture  2 

P = -9,  Ui  = U2  = Q,  I2  = 400-14 

H = Q,  Z = 40.9-14. 

» 

.00 

.021 

.050 

.031 

.046 

.050 

.35 

.037 

.151 

.078 

.102 

.125 

.70 

.096 

.427 

.231 

.336 

.317 

1.05 

.166 

.744 

.423 

.668 

.589 
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Table  4.9  - - continued. 


Amount 

of 

Shift 

Statistics 

F 

RT 

T2 

vn 

nWn 

Mixture  3 

p = .5,  m = [1,  1,  1,  1]T  = -ji2,  = 14,  Z2  = 40014, 

[201.5  1 1 l-i 

o I-  1 201.5  1 i 

^ y’  - 1 1 201.5  1 • 

L 1 1 1 201.5. 

.00 

.035 

.050 

.024 

.033 

.038 

1.40 

.055 

.229 

.047 

.332 

.149 

2.80 

.090 

.444 

.088 

.654 

.216 

4.20 

.172 

.543 

.175 

.769 

.212 

Mixture  4 

P = -9,  Ui  =U2=Q,  hi  =EEt,Z2=  100EEt, 
14  = 0,  E = 10.9-E  Et,  EEt  is  defined  in  (4.1.2). 

.00 

.033 

.072 

.030 

.035 

.053 

1.60 

.076 

.144 

.219 

.292 

.295 

3.20 

.224 

.402 

.617 

.834 

.723 

4.80 

.395 

.738 

.808 

.986 

.880 

Entries  . Proportion  of  times  each  test  statistic  exceeded  the  upper  cc-percentile  of  its 
null  distribution. 

Dimension  : p = 4. 

Sample  Size  : n = 20. 

Number  of  Samples  : rep  = 1000. 
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Table  4.9  - - continued. 

Significance  Level : a = 0.5. 

Distributions  : Quadrivarite  normal  mixture,  choosing  Nj  with  probability  p and  N2  with 
probability  1-p,  where  N;,  i =1,  2,  is  quadrivariate  normal  with  mean  jij  and 
variance-covariance  matrix  Zj.  The  resulting  mixture  has  mean  ji  and 
variance-covariance  matrix  Z (see  (4.1.5)). 


CHAPTER  5 

A MULTIVARIATE  SIGNED  SUM  TEST  BASED  ON  INTERDIRECTIONS 
FOR  THE  ONE-SAMPLE  LOCATION  PROBLEM 


5A  Definition  of  the  Test  Statistic 


We  now  leave  the  repeated  measures  problem,  and  consider  the  one-sample 
multivariate  location  problem.  In  this  section  we  propose  a multivariate  signed  sum  test 
based  on  interdirections,  described  in  section  2.1,  for  the  one-sample  multivariate  location 


Here  0 is  used  without  loss  of  generality,  since  Hq  : £ = fio  can  be  tested  by  subtracting 
fio  from  each  observation  Xj,  and  testing  whether  the  differences  (Xj  - fio)’s  are  located  at 

Q. 

The  procedure  is  somewhat  like  applying  the  multivariate  sign  test,  proposed  by 
Randles  (1989)  for  the  one- sample  location  problem,  to  the  sums 


problem.  To  do  so,  we  let  Xi Xn,  where  X*  = (Xu Xip)T,  be  i.i.d.  as  X = (Xh 

Xp)  , where  X is  from  a p-variate  absolutely  continuous  population  with  location 
parameter  Q (pxl).  We  would  like  to  test 


Ho  : fl*  = Q versus  Ha.-  £*  * 0. 


(5.1.1) 


Xs  + Xt,  1 < s < t < n. 


(5.1.2) 


and  rejects  Hq  for  large  values  of  the  statistic 


(5.1.3) 


where 
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if  (i,  k)*(i’,  k') 


0 


if(i,k)  = (i\  k1). 


(5.1.4) 


and  Cik,i’k'  denotes  the  number  of  hyperplanes  formed  by  the  origin  0 and  p-1  of  the  other 
observations  Xj  (excluding  Xi,  Xk»  Xi’>  and  Xk')  such  that  Xi  + Xk  and  Xi'  + Xk'  are  on 
opposite  sides  of  the  hyperplane  formed.  The  counts  {Qk  i-k-  I 1 < i < k < n,  1 < i*  < k’  < 
n},  called  interdirections,  arc  used  via  7tpik  to  measure  the  angular  distance  between  Xi 
+ Xk  and  Xi'  + Xk1  relative  to  the  positions  of  the  other  observations.  Here  pik  lk.  is  the 
observed  fraction  of  times  that  Xi  + Xk  and  Xi'  + Xk-  fall  on  opposite  sides  of  the 
hyperplanes  formed  by  Q and  other  p-1  observations. 

We  now  examine  some  characteristics  of  the  test  based  on  SS.  First  of  all,  it  can  be 
easily  seen  that  the  test  SS  defined  in  (5.1.3)  is  like  the  multivariate  sign  test  Vn  applied  to 
the  sums  Xs  + Xt>  1 ^ s < t < n.  Secondly,  since  the  p^^.,  1 < i < k < n,  1 < i’  < k’  < n, 

are  invariant  with  respect  to  the  nonsinglar  linear  transformations,  as  shown  by  Randles 
(1989),  it  is  clear  that  SS  is  also  affine-invariant.  Thirdly,  for  p > 1,  SS  does  not  have  a 
small  sample  distribution-free  property.  This  is  because  the  joint  distribution  of  the 
directions  of  Xi  + Xk,  1 ^ i ^ k < n,  from  Q depends  on  the  distribution  of  the  distances  of 
2ii, 1 < i < n,  from  the  origin.  However,  its  large  sample  null  distribution  is  convenient,  as 
is  shown  in  the  next  section.  We  end  this  section  by  proving  a Lemma  which  shows  that 
for  p = 1,  the  test  based  on  SS  is  the  two-sided  univariate  Wilcoxon  signed-rank  test. 
Hence  for  p = 1 it  does  have  a small  sample  distribution-free  property,  as  well  as  a 
convenient  large  sample  null  distribution. 


Lemma  5.1.5  When  p = 1,  the  test  based  on  SS  is  the  two-sided  univariate  Wilcoxon 
signed-rank  test. 
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Proof  of  Lemma  5.1.5  Note  that  when  p = 1, 

1 if  (Xj  + Xk)(Xr  + Xk.)  < 0 

Pik,i'k'  = { 2 if  (Xi  + Xk)(Xi’  + Xk’)  = 0 » 
0 if  (Xj  + Xk)(Xj>  + Xk.)  > 0 


and  thus 


-1  if  (Xj  + Xk)(Xj-  + Xk.)  < 0 
cosC^p^  j.j,,)  = | 0 if  (Xj  + Xk)(Xj-  + Xk0  = 0 
1 if  (Xj  + Xk)(Xj.  + Xk.)  > 0 


Hence,  we  can  write 


cos(KPjkiilc.)  = sign  (Xj  + Xk)-sign  (Xj-  + Xk0 


where  sign(x)  = -1  if  x < 0, 0 if  x = 0,  and  1 if  x > 0. 
Therefore  we  have 

SS  = 4 


—^2  2 X COS(7C  Pik  j'k') 


— 7772  £ 2 sign  (xi  + xk)-sign  (xi*  + xk') 

n(n+D  j<k  j <k. 

~~ 2 { ^siSn  (xi  + xk)}2- 

n(n+i)  j^ 


(5.1.6) 


Noting  that  the  Wilcoxon  signed  rank  test  W+  is  equal  to  the  number  of  positive  Walsh 
averages  ((Xj  + Xk)/2,  1 < i < k < n)  (See,  e.g.,  Randles  and  Wolfe  [1979],  p.57  and 

p.83),  thus  we  have 


X sign  (Xj  + Xk)  = W+  - W, 
i<k 
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where  W'  is  equal  to  the  number  of  negative  Walsh  averages  ((Xj  + X^/2, 1 < i < k < n). 

Now  using  the  fact  that  with  probability  one 

W+  . w-  = Jfo+1)  t 


expression  (5.1.6)  is  equivalent  to 


SS  = 


n(n+l)' 


{W+-  W) 


-l  2 


— - — r { 2W+  - 2 

n(n+l)2  2 


16 


n(n+l)“ 


{w+- 


+ n(n+l)|2 


Au+  n(n+l)^2  n(n+l)(2n+l) 

^ 4 > 24 


n(n+l)(2n+l) 

24 


n(n+l)2 

16 


(W+-EHf,(W+))*  2(2n+1, 
VarH„(W+)  3(n+l)  ’ 


since 


EHo(W*)  = and  VarjyvT) 


(See,  e.g.,  Randles  and  Wolfe  [1979],  p.56.)  This  completes  the  proof. 


Note  that,  when  p = 1, 


— > %2  under  Hq  as  n —>  oo. 


This  follows  from  the  fact  that 


70 


w+~  EHq(W+)  d 

V arH0(W+) 


— > N(0,  1)  under  Hg 


(See  Randles  and  Wolfe  [1979],  p.85),  and 


2(2n+l) 

3(n+l) 


4—  as  n — > oo. 


5.2  Some  Intermediate  Results 


In  this  section  we  present  some  basic  and  important  results  which  will  be  useful  in 
the  next  section  for  finding  the  limiting  distribution  of  SS  under  Hq.  Here  we  discuss  some 
properties  involving  the  sum  observations  Xj  + X^,  1 < i < k < n,  when  the  original  sample 
Xj,  1 < i < n,  is  from  the  family  of  elliptically  symmetric  distributions.  The  first  one  is 
about  the  sum  and  difference  of  two  i.i.d.  spherically  symmetric  random  vectors,  a special 
case  of  elliptically  symmetric  distributions. 

Theorem  5.2.1  Let  Xi,  X2  be  i.i.d.  spherically  symmetric  random  vectors.  Then  X4+X2 
and  Xj-X2  are  spherically  symmetric  random  vectors. 

Proof  of  Theorem  5.2. 1 Let  D be  any  pxp  orthogonal  matrix.  It  suffices  to  prove  that 


D(X1±X2)  = X1±X2 


(see,  e.g.,  Muirhead  [1982],  p.  32),  where  — is  read  “has  the  same  distribution  as”.  Note 


that  we  have 


since  Xj,  X2  are  spherically  symmetric. 

Hence,  by  the  independence  of  Xj  and  X2,  we  have 
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Therefore,  see,  e.g.,  Randles  and  Wolfe  (1979) , p.  16, 

g(X1,X2)  = g(DX1,DX2), 


where  g(x,  y)  = x + y. 

Hence 

Xj  +X2  = DX1  + DX2  = D(X1+X2). 

The  difference  is  handled  similarly.  This  completes  the  proof. 

Next  we  consider  the  general  case  of  two  i.i.d.  elliptically  symmetric  random 
vectors. 

Theorem  5.2.2  Let  Xj,  X2  be  i.i.d.  elliptically  symmetric  random  vectors.  Then  Xj+X2 
and  XrX2  are  elliptically  symmetric  random  vectors. 

Proof  of  Theorem  5.2. 1 Since  Xj,  X2  are  i.i.d.  elliptically  symmetric  random  vectors,  we 
can  write 

X j = A X j and  X2  = — —2  ’ 

* * 

where  Xj,  X2  are  i.i.d.  spherically  symmetric  random  vectors  and  A is  a pxp  nonsingular 
matrix.  Thus 

Xi±X2  = AXj±AX*=A(X*±X*  ). 

This  completes  the  proof  since  Xj±  X2  are  spherically  symmetric  random  vectors,  as 
shown  in  Theorem  5.2.1. 


Define 
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(5.2.3) 


where 


We  examine  the  property  of  hQ^,  X2)  in  the  next  theorem. 

Theorem  5,2.4  LetX.i,X2  be  i.i.d.  spherically  symmetric  random  vectors.  Then 
h(Xi,  X2) is  uniformly  distributed  on  the  p-dimentional  unit  sphere  Sp. 

Propf  of  Theorem  5-2.4  Since  Xi,X2  are  i.i.d.  spherically  symmetric  random  vectors, 
Xj+X2  is  also  spherically  symmetric,  as  shown  in  Theorem  5.2.1.  Thus  it  follows  that 


(See,  e.g.,  Muirhead  [1982],  p.  38,  Theorem  1.5.6.)  The  proof  is  complete. 


Note  that  under  the  assumptions  of  Theorem  5.2.4 


EHoto<*'-  = 2 = EX,(EX2tJl«l.  x2)  I X,]). 


Define 


h*(Xi)  = Exk[h(Xi,  Xk)  I XJ,  1 < i < k < n. 


(5.2.5) 


We  are  interested  in  the  behavior  of  h*(Xj),  1 < i < n. 


Theorem  5,2.6  Let  X1?  X2  be  i.i.d.  spherically  symmetric  random  vectors.  Then  ^(Xj)  is 
a spherically  symmetric  random  vector. 

Prpof  of  Theorem  5,2.6  Let  D be  any  pxp  orthogonal  matrix.  It  suffices  to  show  that 
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D h*(X,)  = h*(Xj). 

(See,  e.g.,  Muirhead  [1982],  p.  32.)  Note  that 
Dh*^)  = DExJhQCLX^IXj 

X + X 

= D EX2[T|--;+-22  n I xj  (by  expression  (5.2.3)) 

_ rPX1+DX2  , 

“ bX2L  H X!  + X2  II  -1-1 

r DX,+DX2 

= EX2L  ||  p x i + d x2  II  * Xi J (II D II  = 1 since  D is  orthogonal) 

= Ex2[h(D  X j , D X2)  I Xj  (by  expression  (5.2.3)) 

= Ex0[h(P  Xj,  X2)  I Xj  ( D X2  = X2  by  spherical  symmetry) 

= h*(DXj)  (by  expression  (5.2.7)) 
d * d 

= h (Xi)  (because  D Xj  = Xj  by  spherical  symmetry).  (5.2.7) 
This  completes  the  proof. 

Hence,  for  i.i.d.  spherically  symmetric  random  vectors  Xj  andX2,  it  follows  that 
* 

we  can  express  h (Xj)  as 

h*(X!)  = Ri  Uj, 

* . * 

where  Uj  ~ Uniform(Sp)  is  independent  of  Rj,  a positive  random  variable.  Defining 

X2  = EHotQi^Xj^VQC1)],  (5.2.8) 


we  get 
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Var-Cov^h*^)]  = 


since 

^ =eHokr;  uI)t(rI  ut)] = %0kr1  )2i- 

In  summary  we  have  the  following  result : if  Xj, ....  Xn  are  i.i.d.  spherically  symmetric 
random  vectors,  then  h*(Xi),  h*(Xn),  defined  in  (5.2.5),  are  i.i.d.  spherically 

symmetric  random  vectors  with 

Ej^QfOfc)]  = Q and  Var-CovHo[h*(X1)]  =^-Ip, 

where  x2  is  defined  in  (5.2.8).  Expressing 

Xj  = Rj  Uj  and  h*(Xj)  = R*  U*,  1 < i < n, 

we  now  try  to  relate  R*  to  Rj  and  U*  to  Uj.  This  is  stated  in  the  next  theorem. 

Theorem  5,2-9  Let  X j , X2  be  i.i.d.  spherically  symmetric  random  vectors  with  X, 
= Rilli»  i = 1,  2,  where  H,  ~ Uniform(Sp)  is  independent  of  Rj.  Then  we  can  write 

h*(Xi)  = R*(R1)H1,  (5.2.10) 


U1)(R1  Uj)1]  (since  =2) 

EH„KRi  )2]EHo[ui(IiI)T] 

Iw®’  )2]ip 


where 
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R1/R2+U21 

r1/r2+u2i)2+i-u21 


(5.2.11) 


is  independent  of  Hi,  andH2  = (U21,  ...,  U2p)T. 

Proof  of  Theorem  5.2.9  Let  X = R U and  h*(X)  = R*  U*.  First  we  want  to  show  that 
U*  = U.  Define 


Ho  = (1,0,  ...,0)T  and  Xq  = R Ho- 


Then 


h*(Xo) 


Kq+X2 
II  2Lo+X2  II 


I Xq]  (by  expression  (5.2.5)) 


= ER2.H2{ 


"R+R2U21" 

R2U22 


L R2U2p  j 


N 


(r+r2u21)2+r2u22 


2 +-+kH 


= eR2.U2{ 


■R/R2+U21' 

U22 


L U2p  J 


N 


(R/R2+U21)2+U22+...+U2p 


7 


R/R2+U21" 

U22 


L u2c 


N 


(R/R2+U2i)2+1-U^ 


} 
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(since  £ U ?•=  1) 
j=l 


" eR2.U2i{eU22,...,U2d[ 


'R/R2+U21' 

U22 


r2 


. u2I] } 


u2o  J 


Noting  that,  given  R2  and  U21,  the  expected  value  of  U2j,j  = 2, p,  is  zero  and  thus  we 
write  the  above  expression  as 


t’CSo)  = eR2,U21{ 


R/^2+U21 

0 


L 0 J 


A^(  R^2+U2l)2+1-U21  } 


= eR2,U2i{ 


= R (R)  Uq, 


R/R2+U21 


(R/R2+U2l)2+1-U21 


} 


LoJ 


(5.2.12) 


where  R is  as  defined  in  (5.2.1 1).  Let  D be  any  pxp  orthogonal  matrix.  Now  we  need  to 
show  that 


h (DX0)=R  (R)  DUq. 


Note  that 

h*(DXo) 


= EX2[h(DXo,X2)IJCo] 

, r D Xp+X2  -i 

■ EX2[  II  D Xn+Xo  II  1 
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r D_  X n+D  X,  -|  d 

= Ex2t  II  D Xq+D,  X2  II  1 (smCC  -2  = - -2) 

r 2.  Xn+ D Xo  -1 

= Ex2L  ||  Xq+X2  11  1 ^ — 11  = 1 since  — 1S orth°g°nal) 

r Xp  + Xo  -1 

■ Hoi'll 

= D h*(Xo)  = D R*(R)  Ho  (follows  from  (5.2.12)) 


= R*(R)DUo. 


(5.2.13) 


Using  expressions  (5.2.12)  and  (5.2.13),  we  have  shown  that  whenever  X = R U , 
h*(X)  = R*(R)  U.  The  indenpendence  of  R*(R)  and  U follows  from  the  independence  of  R 
and  U and  the  fact  that  R*(R)  depends  on  R only.  This  completes  the  proof. 


All  the  above  theorems  and  discussions  lead  to  the  following  result.  Let  Xj, ....  Xn 
be  i.i.d.  spherically  symmetric  random  vectors  with  Xj  = Rj  Hj,  1 < i < n,  where  LLj  ~ 
U(Sp)  is  independent  of  Rj.  Then  h*(X.i),  ...,  h*(Xn)  are  i-i-d.  spherically  symmetric 
random  vectors  with  h (X;)  = R (Rj)  Uj,  1 < i < n,  R (Rj)  is  independent  of  Uj, 

2 

Ej^th*^!)]  = 0 and  Var-Cov^h*^)]  Ip> 
where  R*(Rj)  and  T2  are  defined  in  (5.2.1 1)  and  (5.2.8),  respectively. 


5.3  Asymptotic  Null  Distribution  of  SS 

In  this  section  we  establish  the  null  limiting  distribution  of  SS  under  the  class  of 
elliptically  symmetric  distributions,  with  density  function  as  defined  in  (3.2.1).  Since  the 
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test  based  on  SS  is  invariant  with  respect  to  nonsingular  linear  transformations,  it  suffices 
to  consider 


Xi  = RiUi,  1 <i  <n,  (5.3.1) 

where  Ufs  are  i.i.d.  Uniform(Sp)  random  vectors  and  are  independent  of  R,’s,  i.i.d. 

postive  random  vectors.  As  a first  step,  we  now  seek  an  asymptotic  approximation  for  the 
test  statistic  based  on  SS,  which  possesses  the  following  two  properties  under  Hq  : (a)  the 

difference  between  SS  and  the  approximating  statistic  converges  to  zero  in  probability,  and 
(b)  the  limiting  distribution  of  the  approximating  statistic  is  easily  established.  One 
candidate  satisfying  these  considerations  is 

SS*=^X  I cos(aik,iY),  (5.3.2) 

n(n+i) 

where  aik,i'k'  = the  angle  between  X j + Xk  and  Xj-  + Xk'-  The  first  property  of  this 
approximating  statistic  is  stated  and  proved  in  the  next  theorem. 


Theorem  5.3.3  If  the  observations  Xj’s  are  as  defined  in  (5.3.1)  and  Hq  is  true,  then 


SS  - SS*  — V 0 as  n — > oo. 


Proof  of  Theorem  5.3.3  In  proving  this  result  we  will  need  to  show 


EHq[(SS  - SS*  )2]  — ^ Oasn-^oo 


Note  that,  using  expressions  (5.1.3)  and  (5.3.2), 

.2 


E^KSS  - SS*  )2]  = z ESS  Eh  {[cosfnp^O- 

n (n+i)  i<&i'<k's<ts'<t' 


cos(aik,i'k')Hcos(7t  pst,s't')  - cos(ast,S't')] } 


(5.3.4) 
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It  can  be  easily  shown  that  if  (i,  k)  = (i',  k')  or  (s,  t)  = (s',  t’),  then  the  term  is  zero.  Also, 
if  any  one  of  the  pairs  (i,  k),  (i',  k’),  (s,  t)  or  (s',  t')  is  disjoint  from  the  integers  in  the 
others,  then  the  expected  value  is  zero.  The  outline  of  this  proof  is  given  in  Appendix  B. 
An  examination  of  the  terms  in  (5.3.4)  shows  that,  the  number  of  terms  with  no  disjoint 
pair  is  of  the  order  n6.  This  is  explained  as  follows.  Suppose  pair  (i,  k)  is  not  disjoint  from 
the  integers  in  the  other  pairs,  then  there  exists  a pair  with  at  least  one  common  integer  as  in 
(i,  k).  Without  loss  of  generality,  let  us  assume  this  pair  is  (i,  k').  Now  the  pair  (s,  t)  could 
be  disjoint  from  pairs  (i,  k)  or  (i,  k'),  but  must  have  at  least  one  integer  common  with  that 
of  (s',  t').  (Otherwise  pair  (s,  t)  will  be  disjoint.)  Let  us  call  pair  (s’,  t')  as  (s,  t').  Thus  the 
number  of  terms  with  no  disjoint  pair  is 

v v v v i rc(n+  1 ) _ 1 ) n _ 6 

i<ki'<k's<ts’<t' 

Therefore  we  can  bound  (5.3.4)  by 

(constant)  EHo  { [cos(tc  Piki-k.)  -cos(aikti'k*)]*[cos(7C  pst.s’f)  - cos(ast,s’t’)] } 
< (constant)  EHQ{  \k  Pik>ik.  - ocik,i-k*l  • I7C  pSt,s't'  - ast,s't'l } . 

Note  that 

Pik.i'k'  ‘ aik,i'k'l‘ Itt  Pst.s't'  ~ CCst,s't'l  — 2-2  = 4 a.e. 

Furthermore,  the  n p's  are  consisent  estimators  of  their  respective  a's.  Thus,  by  the 
Lebesgue  Dominated  Convergence  Theorem,  see,  e.g.,  Chow  and  Teicher  (1978),  p.  99, 
we  have 

eHq{|71  Pik.i'k'  • aik,i'k'Mrc  Pst.s’t'  - ast.s’t'l } “ > 0 as  n ->  00. 


Hence  we  have  proved 
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EHo[(SS  - SS*  )2]  — > 0 as  n — > oo. 

This  completes  the  proof,  since  convergence  in  quadratic  mean  implies  convergence  in 
probability. 

>|c 

The  next  theorem  proves  the  limiting  distribution  of  SS  under  Hq. 


Theorem  5.3.5  If  the  observations  Xj,  1 < i < n,  are  as  defined  in  (5.3.1)  and  Hq  is  true, 
then 


SS* 
4 x2 


as  n — > oo; 


where  x2  is  defined  (5.2.8). 

Proof  of  Theorem  5.3.5  Note  that 


aik,i'k’  = angle  between  (Xi+Xk)  and  (Xj'+Xk-) 


= angle  between 


Xi+Xk 

II  Xi+Xk  II 


X p+Xk' 

II  Xi'+Xk’  II  ■ 


Hence,  we  can  write 


cos(aik>i'k0  = ( 


Xj+Xk 


£j'+X.k 

X.i'+2Lk 


since 


Xi+Xjc 

II  Xi+2Lk  II 


has  norm  1. 


Now  SS*  in  (5.3.2)  can  be  written  as 


SS*  = 


4P  y v ( — j+X-k  \T  / Kj'+Xk'  N 

n(n+l)2  ^ i(<k'  " -i+~k  " ; ^ " Ki'+Kk’  " ^ 


= 4P 

n(n+l)2 


(Z-Ii 

i<k 


Xj+Xk  T Xj'+Xk' 

xi+xkii2  '\4l"  -i,+Kk’ " 
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_ I V n X4+X.IC  T v n ^ 2Li’+Xk' 

PV  n(n+ 1 ) ^ II  Xi+Xv  II  ; V n(n  + l)  L II  X;.+Xt.  II  > 

' 2 i^k  ~ 1 ~k  2 i'^k'  ~ 1 ~ k 


= p-R  R , 
v — n— n’ 


(5.3.6) 


where 


R_  = 


V n 

n(n+ 1 ) 
2 


i<k 


Xj+Xk 

"fi+ST 


Moreover  let  us  write 


/ V n Xj+Xk 

k n(n  + l)  Z II  Xi+Xv  II 
5 i<k 


V n Xi+Xk 

n(n  + 1 ) L II  Xj+Xv  II  ; 
2 1=k 


Vn  (■ 


n- 1 
n+1 


■)uln  + e 


2Vn 

n+1 


■)  H2n. 


(5.3.7) 


where 


Xj+Xk 

Xj+Xk  II 


and  U2n  =~^ 


n 

I 

i=l 


Xj 

II  2Lj  II  ‘ 


Since,  under  Hq,  h ^ 1 n , 1 < i < n,  are  i.i.d.  Uniform  (Sp)  random  vectors  (as  shown 
in  Theorem  5.2.4),  we  have 

p 

ll2n  ' 0 as  n — > 00. 

Using  the  fact  that  —>  1 and  — > 0 as  n-+  00,  we  see  that,  under  Hq,  expression 

(5.3.7)  implies 


Rn  - Vn  Uln  —>  0 as  n — > 00. 


(5.3.8) 
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Combining  expressions  (5.3.6)  and  (5.3.8),  now  it  suffices  to  find  the  limiting  distribution 
of  Vii  Ujn  under  Hq.  Define 

llln=  \ ■ Xh*(Xi),  (5.3.9) 

1=1 

$ 

where  h is  as  defined  in  (5.2.5). 

It  follows  from  Theorem  3.3.13  (Randles  and  Wolfe  [1979],  p.  82)  that,  under  Hq, 

Vn  (Uln  - Uln  ) 0 as  n ->  oo, 

and  hence  expression  (5.3.8)  implies 

Rn-^  y;„ 

Note  that,  under  Hq, 

VH  u;n  = VH 

♦ 

where  h (Xj)'s  are  i.i.d.  with 

2 

EHo[h*(X,)]  = Q and  Varj^Ih'QC,)]  =— Ip 
(see  section  5.2).  Hence  by  the  multivariate  central  limit  theorem  we  have,  under  Ho, 

Vn  U*n  Np(0,-y-Ip)  as  n — > oo.  (5.3.11) 

Using  expressions  (5.3.10)  and  (5.3.11),  it  follows  from  Slutsky's  Theorem  that, 

E.n  Np(Q,-y-Ip)  as  n — > oo, 


p 

— ^Oasn— »oo.  (5.3.10) 


( ; • X 2k*(Xi)), 
i=l 
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and  hence 


-^•Rn  Np(0,  Ip)  as  n 


— > oo. 


Therefore,  under  Hq,  using  (5.3.6),  we  have 

-^=(#En)T(#K„)  4X*asn 


4x" 


2x 


2x 


— > oo. 


This  completes  the  proof. 


Now  we  are  ready  to  state  the  theorem  for  the  limiting  distribution  of  SS  under  Hq. 


Theorem  5.3.12  If  the  observations  Xj,  1 < i < n,  are  as  defined  in  (5.3.1)  and  Hq  is  true, 
then 

SS  A 2 , 

—3 > XD  as  n -»  oo, 

4xz  p 

where  x2  is  defined  (5.2.8). 

Proof  of  Theorem  5.3.12  With  Theorems  5.3.3  and  5.3.5  established,  the  result  follows 
directly  from  the  Slutsky's  Theorem. 


Note  that  when  p = 1,  x2  = 1/3.  Thus  when  p = 1, 

SS  3 cc  A 2 
~^T~  SS-^Xiasn->°°. 

Finally,  to  perform  the  test  will  require  a consistent  estimator  of  t",  say  The  test  will 
then  be  based  on  the  fact  that  under  Hq, 
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-- --  as  n — » oo. 

4? 

For  consistent  estimator  of  x , we  consider 

^ ~ n(n-l)(n-2)  ^ ^ cos(7X  Pik,ik’)»  (5.3.13) 

i=l  k<k' 
k kVi 


where  Pjk.ik'  is  as  defined  in  (5.1.4).  The  consistency  is  proved  in  the  next  theorem. 


Theorem  5.3.14  Let  x2  be  as  defined  in  (5.2.8).  Then  x2  is  a consistent  estimator  of  x2 
under  H0. 

Proof  of  Theorem  5.3.14  Define 

n(n-l)(n-2)  ^ ^ cos(aik,ik')  • 

i=l  k<k’ 

k kVi 

Using  the  same  arguments  as  in  the  proof  of  theorem  5.3.3,  it  can  seen  that  under  Hq, 

A 0 as  n -4  oo. 

So  we  can  conclude  our  proof  by  showing,  under  Hq, 

9 P. 

xz  - x — > 0 as  n — > oo. 


Note  that 

't2  = EHn[(h*(Xi))Th*(X1)]  (by  expression  (5.2.8)) 


_ r (X1+X2)t(X1+X3) 

= Ef^L  II  X.i+X.2  II  II  X!+X3  ii  J (by  exPressions  (5.2.3)  and  (5.2.8)) 
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r (Xi+Xk)T(Xi+Xk.)  . 

= EH0[  II  Xj+Xk  II  II  X-i+X-ic'  II  I*  1 - 1 - n>  k < k',  k * l,  k'  * i. 

= Ej^tgCXp  Xk,  Xk0] , (5.3.15) 


where 


fiCXi.  Xk)  Xk0  = 


(Xi+xk)T(xi+xko 

II  Xi+Xk  II  II  Xj+Xk'  II 


satisfies  g(Xj,  Xk,  Xk>)  = g(Xj,  Xks  Xk).  Also  we  can  write 


n(n-l)(n-2)  ^ ^ cos(aik,ik’) 

i=l  k<k' 

k kVi 

2 “ (Xi+Xk)T(Xi+Xk.) 

n(n-l)(n-2)  ^ ^ H X.i+Xk  II  II  Xj+Xk'  II 

k kVi 

2 yyy  f (Xi+Xk)T(Xi+Xk.) 
n(n-l)(n-2)  i<k<k'  l II  Xj+Xk  II  II  Xi+Xk-  II 

(Xk+Xi)T(Xk+Xk.)  (Xk.+Xi)T(Xk.+Xk)  1 

+ II  Xk+Xj  II  II  Xk+Xk.  II  + II  Xk.+Xi  II  II  Xk.+Xk  II  J 

J_  yyy  1 / (Xi+Xk)T(Xi+Xk.) 

(jj)  i<k<k-  3 1 II  Xi+Xk  II II  Xi+Xk.  II 


(2£k+Xi)T(2Ck+Xk.)  (Xk.+Xi)T(Xk.+Xk)  \ 

II  2Ck+Xi  II  II  xk+xk.  II  + II  Xk-+Xi  II  II  Xk-+Xk  II  j 


= QjilS  «<*■*■*>• 


where 

«*(Xi.  Xk,  xk.)  = 


l { (Xi+Xk)T(Xi+Xk.)  (Xk+Xi)T(Xk+Xk.) 

3 l II  Xi+Xk  II  II  Xi+Xk.  II  + II  2Ck+X,  II  II  Xk+Xk'  II 
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(Xk'+Xi)T(X),+Xk)  1 

II  Xk.+Xi  II  II  Xk.+Xk  II  J 


(5.3.16) 


{ |(Xi,  xk,  2£k’)  + g(Xk,  Xi,  Xk)+  g(Xk.,  Xi,  &) } 


is  symmetric  in  its  arguments.  Thus,  by  Corollary  3.2.5,  Randles  and  Wolfe  (1979), 
p.  71,  we  have,  under  Hq, 


x2  - x2 


— ^ 0 as  n — > °°. 


5.4  Asymptotic  Distribution  of  SS  under  Contiguous  Alternatives 

In  this  section  we  establish  the  asymptotic  distribution  of  SS  under  a sequence  of 

$ 

alternatives  approaching  the  null  hypothesis  Hq  : 0 = 0 for  a specific  class  of  elliptically 
symmetric  distributions.  As  a first  step,  let  us  assume  Xi,  ...,  Xn  are  i-i-d.  as  X = 
(Xj, ...,  Xp)  , where  X is  elliptically  symmetric  with  a density  function  fx  of  the  form 

fx(i)=Kpl£rwexp(-[(x-fi’)T£-1(i-fi')/C0]v},xs  RP,  (5.4.1) 

$ 

where  Kp,  Cq  are  as  defined  in  (2.3.2),  £ is  the  point  of  symmetry,  and  E is  the  variance- 

covariance  matrix.  Since  both  SS  and  Hotelling’s  T2  are  affine-invariant,  we  can,  without 
loss  of  generality,  assume  that  I = Ip,  the  pxp  identity  matrix.  Thus,  under  H0, 
X], ...,  Xn  are  i.i.d.  with  density  of  the  form 


fx&)  = Kp-exp{-[(xTx)/Co]V}. 


(5.4.2) 
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Under  the  sequence  of  alternatives  considered  in  section  2.3,  Xq,  X^  ^ i-i-d.  with 
density  fx(x-cn'1/2),  where  fx  is  given  in  (5.4.2)  and  c e Rp-{0)  is  arbitrary,  but  fixed. 

Substituting  p+1  for  p in  the  proof  of  Appendix  A,  we  have  shown  that  if  4v  + p > 2,  then 

r cT9fx(x)  -i  2 

0 < I£(fx)  = } [ fx^)  ] fx(x)dx  < o°  for  all  c ^ 0.  (5.4.3) 

Thus,  when  4v  + p > 2,  the  previous  results  in  conjunction  with  LeCam's  Theorems  on 
contiguity  (see  Hajek  and  Sidak  [1967]  pages  212-213)  can  be  used  to  establish  the 
asymptotic  distribution  of  SS  under  the  sequence  of  contiguous  alternatives.  This  is  stated 
in  the  next  theorem. 


Theorem  5.4,4  Let  Xq,  ...,  Xn  be  i.i.d.  from  an  elliptically  symmetric  distribution  with 
density  function  fx  given  in  (5.4.1)  with  X = Ip,  the  pxp  identity  matrix.  If  4v  + p > 2, 

then  under  the  sequence  of  contiguous  alternatives  for  which  X has  density  of  the  form 
fxOi-cn'1^2),  defined  in  (5.4.2),  we  have 


SS 

4t2 


4v2(-^-f  E 

Xp  Vpcf1 


H0 


[R2V-1R*(R) 


]}V£), 


as  n — > °o,  (5.4.5) 


O * 

where  x and  R (R)  are  as  defined  in  (5.2.8)  and  (5.2.1 1),  respectively. 
Proof  of  Theorem  5.4,4  Define 


1 " T £T9fx(Xi)  ] 

fX^>)  J> 


n * 

Ih  (Xj), 

i=l 


and 
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s„  = vs  iTu;n, 


where  c , X e Rp-  { 0 } , and  h*  is  defined  in  (5.2.5).  Here  X;  = Rj  Uj,  where  Uj  is 

2 T 

distributed  uniformly  on  the  p-dimensional  unit-sphere  independent  of  Rj,  and  R;  = Xj  X;, 


thus  we  can  write 


„ - i fiTSO,T&)^ 

CJVH  i=i 


T = 

1 n — 


2v  i tV'-'a 


CqVH  i=l 


Thus,  under  Hq, 


El- 


— y 


2lTh  (Xj) 
— cTRi2v'1  ijj 

v - i 


— X 

VSfei 


21tR  (Rj)Uj 

-cTRi2v-'  Hi 
C0 


(since  h (Xj)  = R (Rj)  Uj,  as  proved  in  Theorem  5.2.9) 


_J_  " r Ai 
Vn  i=i  L B i _ 


(say), 


where 


Aj  = 2^TR*(Ri)  Ui  and  Bj  = — cTRi2v‘1  Uj. 

rv 

'“O 


Note  that,  under  Hq,  are  i.i.d.  with 
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E 


rAii 

LBiJ 


= Q 


and  Var 


all  °12 

a12  °22-T 


where 

an  = E(A2) 

= E(2kTR*(Rj)  Ui  -2R*(Ri)  U;T2J 
= 4E{[R*(Ri)]}2^TE(UiU7)^ 

= E{  [R*(Rj)]  }2-A,t2-, 
a 22  = E(B2) 

= — ^ E(Ri4v'2)-cTE(UiU7)c 
C0 

= -^-E(Ri4v-2)-cTc, 

PC0 

and 

an  = E(AjBj) 


= E(2X,TR*(Rj)  Uj-  ~ Rj2v*^  Uj7c) 

co 

= ^E(Ri2v-1R*(Ri)).2tTE(UiU7)c 
C0 

= ~~  E[Ri2v'1R*(Ri)]2tTc. 

PC0 


90 


Thus,  under  Hq,  see,  e.g.,  Serfling(1980),  Theorem  B,  p.  28, 


°12 

a22- 


as  n — » oo. 


where  a^,  (722,  and  are  as  defined  above.  Applying  LeCam’s  Theorems  on  contiguity 
(see  Hajek  and  Sidak  [1967],  p.  212-213),  if  4v  + p > 2,  then  we  have,  under  the  sequence 
of  alternatives, 

d 

Sn  — > N(o12,  an)  as  n — > oo. 


That  is,  under  the  sequence  of  alternatives, 


Sn  = Vn  lTu\n  4 n(-^  E[Ri2v'1R*(Ri)]-^T£, 


pC, 


4r 

P 


■lTk) 


as  n — > oo 


since 


E{[R*(Ri)]}2  = T2. 


Therefore,  under  the  sequence  of  alternatives, 

•fr,  ujn  4 NP(^;  E[Rj^v''R*(Ri)]  c,  ^ Ip)  as  n -> 
Pc0 


and  hence 


^VSu;n4Np(741;E[R>-1R*(Ri)].£.Ip)  as  n 
V ptCq 


— > oo. 


Recall  that 


Rn  - Vn  U1n  *4  0 as  n — » oo  and  SS*  = p-RTR 

11  — In  r — n — n 
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(see  expressions  (5.3.8)  and  (5.3.6),  respectively).  Using  the  Slutsky’s  Theorem  we  see 
that  Rn  and  Vn  U ln  have  the  same  limiting  distribution,  and  thus,  under  the  sequence  of 

alternatives,  we  have. 


SS 

4x2 


V*ta 


4 x 


= jt-A) 

2x  ^ 2x  ^ 


“^Xd('T^;{EhJr2V  1r*(r)]  } Vs),  asn 


“p  V xVf  H° 


— * oo. 


Therefore,  under  the  sequence  of  contiguous  alternatives. 


SS 

4t2 


4y2(^-{e 

^ V 9 „2v  l n 


^pN2pcn2v 


H0 


[r2v_1R*(R) 


]}Vs), 


as  n — > oo. 


since 


SS 

4x2 


SS*  An 

— 7 0 as  n — > oo 

4x2 


(proved  in  Theorem  5.3.3).  This  completes  the  proof. 


Finally,  we  establish  the  Pitman  asymptotic  relative  efficiency  of  the  test  based  on 
SS  relative  to  Hotelling’s  T2  in  the  next  theorem. 

Theorem  5-4.6  Let  Xj,  ...,  Xn  be  i.i.d.  from  an  elliptically  symmetric  distribution  with 
density  function  f^  given  in  (5.4.1)  with  Z=  Ip,  the  pxp  identity  matrix.  If  4v  + p > 2, 
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then  the  Pitman  asymptotic  relative  efficiency  of  the  test  based  on  SS  relative  to  Hotelling’s 
T2  is 


ARE(  ^4 , T2  ) = 
4xz 


{eHo[r2v-‘r*(R)]}2, 


(5.4.7) 


where  x and  R (R)  are  as  defined  in  (5.2.8)  and  (5.2.11),  respectively. 


Proof  of  Theorem  5.4,6  Since  under  the  sequence  of  contiguous  alternatives. 


the  result  follows  by  taking  the  ratio  of  the  noncentrality  parameter  of  the  test  SS  to  the 


In  this  section  I describe  the  numerical  evaluation  of  the  ARE  expression  given  in 


(RCI)  account  offered  by  the  Northeast  Regional  Data  Center  (NERDC),  running  on  an 
IBM  computer  using  fortran  77.  A subroutine  called  ADAPT  (see  Genz  and  Malik  [1980]) 
was  used  for  numerical  integration  over  an  N-dimensional  rectangular  region. 

Now  we  briefly  describe  how  the  calculations  were  performed.  First  note  that, 
following  some  standard  transformations  and  arguments,  we  can  express 


noncentrality  parameter  of  the  test  T . 


5.5  Numerical  Evaluation  of  AREfSSMx2  T2) 


Theorem  5.4.6.  All  calculations  were  performed  with  a Research  Computing  Initiative 


Kp7iPco+V"1/2 


1 1 

j j f(x,  y)  dxdy, 
0 0 


(5.5.1) 
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and 


x1 2 


Kp7tp/2Cp/2 

v-F(|) 


v ill 

) III  g(x,  y,  z)  dxdydz, 
0 0 0 


where 


f(x,  y)  = 


(-/nx) 

{-[ny) 


1 p-1  p-2v 

]2v)-  (-[nx)  2v(./ny)  2v 


(5.5.2) 


g(x,  y,  z)  = k([  f ^ ]2v)  k([  ) Lrny  ]2v)-(-/nx)  2v  (.(ny)  2v  (.foz)  2v 


(-/■nx) 


p-2v 


E^y  E^y 


(-^y) 


(-/"nz) 


K s+cos(8) sinp  (8)T(^-) 

0 (s2+2s-cos(8)+1)1/2  r(^)7i1/2 


(5.5.3) 


Kp  and  Cq  are  as  defined  in  (2.3.2),  and  In  is  the  natural  logarithm.  Using  expressions 
(5.5.1)  and  (5.5.2),  we  can  simplify  expression  (5.4.7)  to 

[ | j f(x,  y)dxdy]2 

ARE(SS/4t2,  T2)  = constant— j-y-j , (5.5.4) 

j j j g(x,  y,  z)dxdydz 
0 0 0 


where 

4v2n^) 

constant  = — — — — , (5.5.5) 

p2r2(^) 


and  f,  g are  as  defined  in  (5.5.3).  Now  using  subroutine  ADAPT  with  a relative  error 

1 1 

ESTREL,  we  can  approximate  j j f(x,  y)dxdy  by  VAL  satisfying 

0 0 

1 1 

I j j f(x,  y)  dxdy  - VAL  I < ERROR, 

0 0 
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where  ERROR  = VAL  x ESTREL.  Hence  we  approximate  [ j j f(x,  y)  dxdy]2  by  (VAL)2 

0 0 

with  the  maximum  error 

ERROR  1 = 2VALERROR  + (ERROR)2.  (5.5.6) 

1 1 1 

Also,  we  approximate  j j j g(x,  y,  z)  dxdydz  by  TAU  satisfying 
0 0 0 

1 1 1 

I j j j g(x,  y,  z)  dxdydz  - TAU  I < ERROR2,  (5.5.7) 

0 0 0 


where  ERROR2  = TAUxESTREL.  So  we  can  write  expression  (5.5.4)  as 


ARE(SS/4t2,  T2)  = constant-^ 


(VAL)2  ± ERROR  1 , 
TAU  ± ERROR2  ’ 


and  hence  approximate  ARE(SS/4t2,  T2)  by  constant- ((VAL)2/TAU)  with  the  maximum 
error 


ERREST  = constant 


ERROR1TAU  + ERROR2(VAL)2 
TAU(TAU  - ERROR2) 


(5.5.8) 


Note  that  when  p = 1,  using  the  expression  (5.4.6)  on  p.166  of  Randles  and  Wolfe  (1979), 
the  ARE  expression  in  (5.4.7)  can  be  simplified  to 


ARE(SS/4x2,  T2)  = 


12v2T(^) 

21/vr3(i) 


(5.5.9) 


In  Table  5.1  we  present  the  asymptotic  relative  efficiency  of  the  interdirection 
signed  sum  test  relative  to  Hotelling’s  T for  selected  values  of  v and  p.  When  p = 1,  the 
ARE’s  are  evaluated  from  expression  (5.5.9).  For  p = 2 to  5,  and  thus  4v+p  > 2 for  all  v > 
0,  the  ARE’s  are  numerically  evaluated  using  expressions  (5.5.4)  to  (5.5.7).  When  the 
underlying  population  is  multivariate  normal  (v  = 1.0),  Hotelling’s  T2  performs  better,  yet 
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the  test  based  on  SS  appears  to  be  quite  competitive.  For  lighted-tailed  distributions  (v  = 
2.0  to  5.0),  Hotelling’s  T2  still  has  better  power.  For  heavy-tailed  distributions  (v  = .75, 
.50,  .25,  and  .20),  the  interdirection  signed  sum  test  is  more  effective  than  Hotelling’s  T . 
In  Table  5.2  we  display  the  values  of  ERREST,  defined  in  (5.5.8),  which  bound  the  error 
in  the  estimates  of  ARE(SS/4x2,  T2)  for  p = 2 to  5. 
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APPENDIX  A 


In  the  context  of  theorem  2.3.10,  we  show  that  if  4v  + p > 3,  then 

r cT3fv(z)  -.9 

0 < I£(fz)  = J [ ' J fz(z)dz  < 00  for  all  c * 0. 


Proof  Under  Hq  : Q = Q,  taking  A £ AT  = Ip.i,  Z has  density  function 

fz(z)  = Kp  f exp  { - [ ^ q+S  -] V } ds, 


and  thus,  following  standard  arguments,  we  have 

2fzfe)  = Kp  @ eXp{-[-%i]v}ds. 


Therefore, 

yfz) 


rcT3fz(Z)  -I, 

L fz(Z)  J 


rT'y  , „2 


fz(Z)^ L Co 


0) 


Now,  it  suffices  to  show  that 


0 <E(  h§|)]2 [(j°(-a^i)v-1  exp(.e^i)v}  as]2}  < 


T-7 . 02 


(A.l) 


( I ) We  first  show  the  above  expected  value  is  finite.  This  is  done  by  considering  the 
following  two  cases : 
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Case  1 : If  v < 1,  then 


E ( -p(-^fi)v) d^2} 

s E ( [-S)]2  r4r]2v'2[[ ^ 

= %M^2^)2v~2} 

= [^]2c202ve{£tuut£}e{r4v'2} 


T7_lc2 


(This  step  follows  from  the  fact  that  Z = RIJ,  ZTZ  = R2,  and  the 


independence  of  R and  1J.) 


[^-]2c20  2''e(^A}e{r4v-2}. 


■P-1 


So  now  it  suffices  to  show  that  E[R4v'2]  < Note  that  R2  has  density  function  of  the 
form 

Kn7t(p‘1)/2 


fR2(D  =- 


F[(p-l)/2] 


. r(p-l)/2-l  | exp  { -[(r+s2)/C0J V } ds. 


Thus 


e[r4v'2]  = e[(r2)2v'1] 


9V  ir(P-D/2  oo  oo 

Tkp-D/Z)  JJ  1+(P- 1)/2*  1 ,exp { - [(r+s2)/Co] v ) ds  dr 
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9 

Letting  r = t , we  have 

{„  -i  ak 

r4v-21  ( [ ,4v+p-2 exp{-[(t2+s2)/C0]V)ds dt. 

r[<p-i)/2]  (jj 


Letting 

s = VCo  x1/,vsin(0)  and  t = ^fCQ  xlyVcos(0), 

we  have  [(s2+t2)/Q)]v  = x2  and  the  jacobian  of  the  transformation  is  - Cq  x^  V^/V. 

v 

Thus,  the  double  integral  is  equal  to 


n/2 


i i 


[CQ1/2x1/v  cos(0)]4v’4+p  { (l/v)C0x(2'v)/vexp(-x2)}  d0dx 


.71/2 


= (q5vtp/2'Vv  ) { j (cos(e))4v+p'4de}  { j x<3v+p':2>/v  exp(-x2)dx  } . 


i 


Note  that  the  first  integral  in  the  above  expression 


71/2 


J 


(cos(0))4v+p'4  d0  < oo. 


if  4v+p-4  > -1,  i.e.,  if  4v+p  > 3.  Also,  letting  y = x“,  we  have  the  second  integral 
"x(3v+P-2)/v  eXp(-x2)  dx  = I Ty(4v-2+p)/2v-l  exp(.y)  dy  < 00> 


i 


J 


if  (4v-2+p)/2v  > 0.  Note  that  4v-2+p  is  always  positive  since  we  consider  v > 0 and  p > 2. 
So  we  have  proved  that,  if  0 < v < 1 and  4v+p  > 3, 

E { R4v'2  } < oo,  (A. 2) 

and  hence  I£(fz)  < oo. 
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That  is 


if  p > 3 and  0 < v < 1 
if  p = 2 and  ^ < v < 1 


(A. 3) 


Case  2 : If  v > 1, 

then  we  can  write  the  integral  in  (A.l),  provided  it  exists,  as 
eXp ds 


1 7T7 , 2 

co 


ZTZ+s2  \vi  , °°/  ZTZ+s“  v- 1 r / ZTZ+s2  \vi  j 

) } ds+|  ( (^ — V exp{-(-go — )v]  ds 


= A + B, 


where 


A = 


1 ZTZ+s2  w.i 


ZTZ+s2  W1 


J(_Q,  )V'  exPf-(“c0  > >ds 


^ ^ i ZTZ+1  \v-i  r f ZTZ+s“  \vi  j / • 1, 

< J (-^  ^ — ) exp  { -(  ~ ) } ds  ( since  v > 1) 

t Z.^Z.+ 1 w—  1 ^ r / ZTZ+s2  1 j 

= (— cT_)  j expl~(  Co  ) ids’ 


and 


B 


°°  ZTZ+s2  NV.i 


ZTZ+s2  w 


- / « P{-e^^)v)ds 

S j s(-^T|+s2  )v'1  exp{-(  )v}  ds  ( since  s > 1) 
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C0  °°,  ZTZ+s2  w-l  , 2s  x f / ZTZ+s2  wi 

T / (— cT")  < q > exP*-<  c0  > 1 ds 


Co  r f / Z Z+S  i l°°i 
-[-exp(-(-  c— ) }lj  ] 


co  f/ZTz+iw1 

— exp{-( — ^ — ) }. 


2v 


Co 


Note  that 


zTz+s2  1y1 


fz(Z)  = 2KP  J exP(~[  ~ Cq  a ] )ds 
> 2Kp  j exp{-[  ~ ] v } ds 


1 ZTZ+1  -J  V 1 

= 2Kpexp{-[-r7J1  3V). 


> 2Kpjexp{-t^Li-]v}ds 
zTz+l 


C0 


Therefore 


A < |-  zTz+l  jy-i  1 B .C0  1 


fz(Z)  ~ L Cq 


and  r (ry\  — ’ 

2Kn  fZ©  2v  2Kr 


Thus,  by  (A.l),  we  have 


E{t^i)]2[A+B]2} 


^E{[(£Ta2[ti^tL)v-.+^]2} 


^-rE{£TiJUT£}E{R2[e 

4Kn2  1 J L 


Rz+1 


C0 


-)V"1+—  l2  | 
2v  J 
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(Since  Z = RU,  ZTZ  = R2,  and  R is  independent  of  and  U.) 

So  now  it  suffice  to  show  that 

e{r2[<^)v-‘+^]2}<~. 

Let 

g(R2)  . R2[e^L)V-,+  So  ]2  . { R[f^±L)v-.+  So  ] } 2. 

M)  2v  M)  2v 

Note  that  g is  a function  of  R of  order  2v  -1.  Using  the  result  that 


“ If  X e £p,  then  X e for  all  0 < q < p.”, 

it  suffice  to  show  that  e[(R^)^v'^]  < From  (A. 2),  we  see  that  the  result  holds  if 
4v+p  > 3.  Note  that  4v+p  > 3 is  always  true  since  we  are  now  considering  v > 1 and  p > 2. 
Thus, 

Ij.(fz)<~  if  v > 1.  (A.4) 


Combining  (A.3)  and  (A.4),  we  have  proved  I£(fz)  < 00  if  4v+p  > 3.  That  is 


i£(fz)  < ~ 


if  p > 3 and  v > 0 
if  p = 2 and  v > ^ 


(A. 5) 


( II ) We  now  need  to  show  I£(fz)  > 0. 
From  (A.l),  it  suffices  to  show  that 


it  / T cTZ  12  r °°/  ZTZ+s2  w-i  r / ZTZ+s2  \vi  j 12  1 a 
E i LWJ2  [ J (~Ca  ) exP{-(~^  )}  ds]2  ) > o. 


(3 


C0 


Co 
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00  7^74-s^  7T7  , „2 

Since  fz(Z)  > 0 and  h(Z)  = j (-~  )v_1  exp{-(  ^ )v}  ds  > 0 V Z,  we  need 


to 


show 


E{ (cTZ)  g(Z) }2  > 0 V c * 0,  where  g(Z)  = h(Z)/fz(Z).  (A.6) 

To  verify  (A.6),  we  note  that  if  3 c * 0 s.t.  E{(cTZ)  g(Z)}2  < 0 

^ {(£TZ)  g(Z))2  = 0 a.e. 

=7  (£TZ)  g(Z)  = 0 a.e. 

=>  (£TZ)  = 0 a.e.  (since  g(Z)  > 0 V Z) 

=>  (£TLD  = 0 a.e.  (since  Z = RU  and  R > 0) 

=>  Var(cTU)  = 0 

=*  c?  Var(Uj)  = 0 
j=l  J 


(Since  cT  = (cj,  ...,  Cp^),  UT  = (Ult  ...,  Up.!),  and  Uj's 
are  uncorrelated.) 

=>  Cj  = 0 V j = 1, ...,  p-1.  A contradiction. 


This  completes  the  proof. 


APPENDIX  B 


Here,  we  prove  in  the  context  of  Theorem  5.4.4  that  if  any  one  of  the  pairs  (i,  k),  (i’,  k'), 
(s,  t)  or  (s',  t')  is  disjoint  from  the  integers  in  the  others,  then  the  expected  value 

E^UcosOt  Pik,i'k') ' cosCocik.i’k')]  * [cos(tc  pst.s't’)  - cos(ast,s't')] } 


is  zero. 

Proof  Suppose  (i,  k)  is  unique  (i.e.  disjoint),  then 

E{  [cos(7tpik,i-kO-cos(aik>i-k')]-[cos(7tpstjS't')-cos(asts-f)] } 


E ( EXi,  Xkftcos(7tPik.i'k')-cos(aik,i  v)]  ’ 


[cos(7tpst  s't')-cos(ast>s’t0]  I all  Xj  except  j = i and  k)  } 


(B.l) 


Note  that 

EXj,  xj  [cos(^ik,iv)-cos(aikti.k.)]- 
[cos(7tpst>s't')-cos(ast)S't')]l  all  Xj  except  j = i and  k} 

= EDj,  Aj,  Dk,  Ak  t [cosCTtpCDjAi,  DkAk,  XiS  Xk-))-cos(a(DiAi,  DkAk,  Xjs 


Xk’))]-[cos(7tpstyS't')-cos(astiS't0]l  all  Xj  except  j = i and  k} 


= EAj,  Ak  { EDi,  Dk[c°s(^p(DjAi.  DkAk,  Xr,  Xk’))-cos(a(DjAi,  DkAk,  Xr, 
Xk'))l  Ai,  Ak]-[cos(jtpst>s-t')-cos(asts't')]l  all  Xj  except  j = i and  k}  ,(B.2) 
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where  Dj  = signCX^),  Aj  = sign(Xil)Xi.  Now  conditioning  on  Aj,  Ak  and  all  X,  with  j * i 
or  k,  we  have 

EDi,  Dk[C0S(a(DjAi,  DkAk,  Xr,  Xk'))] 

= { cos(a(Ai,  Ak,  XiS  Xk'))  P(Dj  = 1,  Dj,  = 1) 

+ cos(a(Aj,  -Ak,  XiS  2CkO)  P(Di  = 1,  Dk  = -1) 

+ cos(a(-Ai,  Ak,  Xj-,  XkO)  P(Di  = -1,  Dk  = 1) 

+ cos(a(-Ai,  -Ak,  Xr,  Xk.))  P(Dj  = -1,  Dk  = -1) } . (B.3) 

Note  that  D,  and  Dk  are  conditionally  i.i.d.  with  P(D,  = 1)  = P(Dj  = -1)  =^.  Note  also  that 

cos(a(Aj,  Ak,  XiS  XkO)  = -cos(a(-Ai,  Ak,  XiS  XkO), 

and  cos(a(-Ai,  -Ak,  Xi',  Xk-))  = -cos(a(Ai,  -Ak,  Xi',  Xk'))- 

Thus,  the  expression  in  (B.3)  equals  0.  Similarly, 

cos(7tp(Aj,  Ak,  Xis  Xk-))  = -cos(7tp(-Ai,  Ak,  Xi',  Xk0), 

and  cos(jtp(-Ai,  -Ak,  Xi-,  Xk0)  = -cos(7tp(Ai,  -Ak,  Xr,  Xk0), 

since,  for  example,  p(-Ai,  Ak,  Xi',  Xk') = 1-pCAi,  Ak,  Xi',  Xk0.  Therefore,  following  from 
the  same  arguments  made  above,  conditioning  on  Aj,  Ak  and  all  Xj  with  j * i or  k, 

EDj,  Dk[cos(icp(DiAi,  DkAk,  Xi-,  2Ck  ))]  = 0.  (B.4) 

Hence,  using  expressions  (B.2)  - (B.4),  expression  (B.l)  implies  that 

eHq{  Pik,i'k') ' cos(aik,i,k,)]-[cos(7t  Pst.s't')  - cos(ast,S't')] } = 0. 


This  completes  the  proof. 
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